SOCIETY 


ENGINEER 


ERICAN 
CIVIL 


OF 


; 


1981 


4 of Civil Engineers 


0088-8001 
ore 


107 NO. ST4. APRIL 


= 
avil 
BEE | 
OCIETY 
q aa 


OF CIVIL ENGINEERS 


Late 


BOARD OF DIRECTION 
President 


Ievan F. 


President 
Joseph S. Ward 


President Elect 


Vice Presidents a 


Robert D. 
Francis J. Connell 


d Directors 
G. Abegg 
Flo d A. Bishop 

ary Byrd 
tae Feeser 
John A. Focht, Jr. 
Sergio Gonzalez-Karg 


Lyman R. Gillis 


Paul R. Munger 
William R. Neuman 
Leonard S. Oberman 
John D. Parkhurst 
Celestino R. Pennoni 
Robert B. Rhode 
James E. Humphrey, Jr. 2," S. Russell Stearns _ 
Stafford E. Thorntor 


Richard W. Karn William H. Taylor 
Robert E. 


Leon D. Luck 
Eugene Zwoyer, Executive 
Director 
Louis L. Meier, a Counsel/ Assistant 
William H. Wisely, Executive Director Emeritus 


Michael N. Saigo, Treasurer 
Elmer B. Isaak, Treasurer 


» 


Donald A. Buzzell, Managing Director for 
Education and Professional Affairs ¢ 
Robert A. Crist, Jr., Managin Director for 
= Publications and Technical Affairs 
Alexander Korwek, Managing Director for 
Finance and Administrative Services 
Alexandra Bellow, Director, Human 
David Dresia, Director, Publications Poe 
Production and Marketing a 
Barker D. Herr, Director, Membership 
Richard A. Jeffers, Controller 
Cari E. Nelson, Director, Field Services _ 
Don P. Reynolds, Director, Policy, 


and Public Affairs 


David Dresia, Director, Publications 


James M. Director, Public 


Albert W. Turchick, Director, Technical 


Communications 
Services 


George K. Wadlin, Director, Education 
Services 


R. Lawrence Whi ple, Director ee 
Management Services 


Stafford E. Thornton, 
Richard W. Karn 


= G. Abeg 
John A. Focht, Paul R ‘Munger 
William R.Neuman 


STRUCTURAL DIVISION 


Executive 


John E. Bower, Chairman — 

Ronald G. Domer, Vice Chairman 

Peter B. Cooper . 

Donald McDonald, Secretary 

Frederick H. Sterbenz, Management 
Contact Member 

Joseph H. Appleton, Group 
Contact Member 


Gow 


John E. Bower, Chairman and Exec. Comm. 
Julie E. Gibouleau, Assistant to the ana 


Contact Member 
Mihran S. Agbabian 
Marvin E. Criswell 
John T. DeWolf 
Douglas Foutch 
Ovadia E. Lev 


ote Le-Wu Lu 
Bruce E. Lyons 
Walter Podoiny 
Simiu 
James T. P. Yao 


SERVICES DEPARTMENT 


Production and Marketing 


Technical and Professional Publications _ 


Richard R. Torrens, Mana ec 
Joseph P. Cerami, Chief Copy Editor 
Linda Ellington, Copy Editor 
Thea C. Feldman, Copy Editor 
Meryl Mandle, Copy Editor __ 
Joshua R. Spieler, Copy Editor 
Shiela Menaker, Production Co-ordin 
Richard C. Scheblein, Draftsman 


d 
| 
4 

4 
ea 

q Information Services 4. 

"Bruce ickerson, Director, Legisiative services Elan Garonzik, Editor 


| PERMISSION TO PHOTOCOPY | JOURNAL PAPERS 


Permission to photocopy for or internal ‘beyond 
in Sections 107 and 108 of the U.S. Copyright Law is granted by the American | 
_ Society of Civil Engineers for libraries and other users registered with the Copyright ‘ 
_ Clearance Center, 21 Congress Street, Salem, Mass. 01970, provided the appropriate . 
fee is paid to the CCC for all articles bearing the CCC code. Requests for special — 


comes 4 


Nonlinear Finite Element Analysis of P Pull-Out Test = 
Elastic- Plastic Large Displacement Analysis of Pipes . 
_ by A. F. Saleeb and W. F. Chen — 
Intermediate Stiffeners for Thin-Walled Members 


‘ 


by Dann H. Hall . 


od 


_ This Journal is published monthly by the American Society of Civil Engineers. Publications : 
office is at 345 East 47th Street, New York, N.Y. 10017. Address all ASCE correspondence 
to the Editorial and General Offices at 345 East 47th Street, New York, N.Y. 10017. 
Allow six weeks for change of address to become effective. Subscription price to members 
is $22.50. Nonmember subscriptions available; prices obtainable on request. Second-class — 
postage paid at New York, N.Y. and at additional mailing offices. ST. 2 sits 
Ts Society is not ew for any statement made or opinion — in tad 


permission OF u copying shou © a ressed to the anager echnical an 
— 
} 
&g 
— 
— 
689 
| 


Wind- -Excited of Buildings in Higher Modest 


Effect in Slabs: Rigid-Plastic by Mikael 


of Structural il Models,” by C Charles A. Miller 
by Phillip L. Gould. . . 
Weibull Distributions and Extreme Wind by Emil Simiu 


represent part of ASCE | in the EJC | 
: od plan. The retrieval data are placed herein so that each can be 
cut out, placed on a 3 x 5 card and given an accession number for the 
user’s file. The accession number is then entered on key word cons 
so that the user can subsequently match key words to choose the ie 
he wishes. Details of this program were given in an August, 1962 article 


in CIVIL ENGINEERING, | reprints of which are available on request — 


a “Discussion period closed for this paper. Any other discussion received during this _ 
discussion period will be published in in 5 Journals. 


Prog, Paper — 
Limit Design of Extended End-Plate Connections, by Allan P.Mann 
and Linden I Morris (Mar 1979 Prior Discussion: Feb 1980) 
| 
4 
: 
4 
| a Chi st 


16197 ANALYSIS OF /PULL- L-OUT T TEST a 


Concretes; Cracking; Failure (mechanics); Finite element method; Materials sent” 


| om Mathematical models; Matrix methods; Mechanical properties; Pull- 


| out tests; Softening; Strain hardening; Strength; Stresses; Structural 

| analysis; Triaxial stresses wt 

ABSTRACT: A specific pull-out test used to determine in- concrete 

‘Strength is analyzed. This test consists of a steel disc that is extracted from the 

' structure. The finite element analysis considers cracking as well as strain hardening 

; | and softening in the pre- and post-failure region, respectively. The aim is to attain a _ 
| Clear insight into structural behavior. Special attention is given to the failure mode. " 
| Severe cracking occurs and the stress distribution is very inhomogeneous. However, _ 

_ | large compressive forces run from the disc in a rather narrow band towards the support 

and this constitutes the load-carrying mechanism. The failure is caused by the crushing 

of the concrete in this region, and not by cracking, 

| REFERENCE: Ottosen, Niels Saabye, “Nonlinear Finite Element Analysis of Pull- Out 
Test,” Journal of the Structural Division, _ ASCE, =~ 107, No. ST4, Proc. Paper 


Finite element method; Friction; Offshore structures; Pipes; Stiffness; 
ABSTRACT: A numerical procedure based on the finite segment method has been 
presented and applied to large displacement analysis of pipe wnethod bes been 
: eee encountered in offshore structural engineering applications. In this method, — 
_ the pipe is treated as an assemblage of segments and the moment-curvature relation of _ 
_ the tubular cross section is approximated by simple mathematical expressions. Material 
and geometrical nonlinearities are treated using a modified incremental tangent 
_ stiffness approach. A number of numerical check problems are solved and compared 
with available “exact” solutions. Finally, the method is applied to an offshore pipe 


_ REFERENCE: Saleeb, A. F., and Chen, W. “Elastic-Plastic Large Displacement 


Analysis of Pipes,” Journal of the Structural Division, ASCE, Vol. 107, No. ST4, Proc. B: 


WALLED MEMBERS 


5 KEY WORDS: Beams; Cold-rolled steel; ‘Dee Critical point; Plates” 
(structural members); Rigidity; Stability; Stiffeners; Stiffness; Strength of 
ABSTRACT: Am effective width approach is presented for 
-" strengths of thin-walled compression elements with intermediate stiffeners. The study 
_ is limited to intermediately stiffened assemblies in which the stiffener is located in the 
center of the stiffened plate and is of relatively small width compared to the plate. In 
' ry conjunction with procedures for predicting ultimate strength, > a requirement for 
adequate stiffener rigidity is given. Stiffener adequacy is assessed as that rigidity for 
' - which the ultimate strength of each component plate element, one on either side of the 


| stiffener, equals that of an element of similar dimensions and material properties, but — 


‘supported by webs at both longitudinal edges. Critical and post-critical behavior are 

studied analytically and experimentally; however, the stiffener requirement 
' _ developed based primarily on the experimental results. Simple formulations for — 
assessing the performance the stiffener and plate elements are presented. 


é 
REFERENCE: ‘Desmond, Thomas P., , Pekoz, Teoman, and Winter, George, 
| “Intermediate Stiffeners for Thin-Walled Members,” Journal of the Structural Division, 
ASCE, Vol. 107, No. Proc. Paper 16194, 1981, PP. 627- 648 
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16198 PROPOSED STEEL COLUMN STRENGT ITERIA 

KEY WORDS: Columns (supports); Elasticity; Equations; Load factors; 

Plasticity: Regression analysis; Steel; Strength; Strength of materials; 

Ss 


ABSTRACT: The dev elopment of proposed design aaa is presented for pinned end 
columns for the load and resistance factor design provisions being developed for the © 
American Institute of Steel Construction (AISC). Available test data have been 
examined by regression analysis to determine a practical column strength equation for 
design. The study resulted in the division of slenderness into three regions: a plastic 
region where slenderness does not effect column strength; an inelastic region where — 
strength decreases linearly with increasing slenderness; and an elastic region where | 
strength is defined by a hyperbola. Equations representing column strength in each of : 
these regions are developed, and a resistance factor to be used with these equations: is 
7 dy REFERENCE: Hall, Dann H., “Proposed Steel Column Strength Criteria,” Journal of ‘4 
the Structural Division, ASCE, Vol. 107, No. ST4, Proc. Paper 16198, April, 1981, 


16205 EGYPTIAN GYPSUM FOR HOUSING 
KEY WORDS: _ Composition; 


Cairo University the Massachusetts Institute of under the 
Adaptation Program. This program is supported by the Agency for International 
Beg of the United States Department of State. Newly-found large deposits 

- make gypsum attractive. Strength, density, retarders, optimum composition, and — 

_ B i additives such as sawdust and foaming agents are examined. Solid and cored gypsum > 

_ block, and the strength characteristics of walls made of such units, are compared with © rt. 


14 


REFERENCE: Reda Youssef, Mahmoud A., and Dietz, Albert G.H., “Potentialities of 
| Esyptian Gypsum for Housing,” Journal of the Structural Division, ASCE, Vol. 107, 
ST4, Proc. Paper 16205, April, 1981, pp. 671- 
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16211 LATERAL- BUCKLING OF TAPERED I- BEAMS 
KEY WORDS: Buckling; Computation; Energy methods; Finite difference 
theory; I beams; Matrices (mathematics); Stability; Stability analysis — po 
ABSTRACT: The lateral-torsional buckling behavior of non-prismatic beams is not — 
_ fully documented. Using the total change in potential energy during buckling, a fourth- 
_ order differential equation in the angle of twist, B, is obtained, in which the location of 
the transverse load with respect to the centroid is considered. Using central finite 
differences, this equation is transformed to a characteristic matrix eigenvalue equation. 
_ Solutions for a number of conventional tapered-beam problems are presented. — 
_ REFERENCE: Brown, Thomas G., “Lateral-Torsional Buckling of Tapered I-Beams,” 1? 
_ Journal of the Structural Division, ASCE, Vol. ae ST4, Proc. Paper et, 


| —_ 
(| 
| 
j 
| 
4 
| 
f 
4 
= 


JOURNAL OF THE 


fod 
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A considerable ir interest is in-situ proper- 
ties, and various destructive as well as nondestructive methods are currently 
applied. Knowledge of the in-situ concrete compressive strength is of particular — aq 
importance and pull-out tests have been proposed quite early for this purpose E 


(14,15). Newer investigations within this field have been reported, €. by 
; Malhotra (9), Clifton (1), and Mailhot, et al. . (8). The importance of those pull- out 
tests, in which a circular steel disc is extracted from the structure using . 

; _ cylindrical counter pressure, has manifested itself as a tentative standard, i.e., 


‘The present paper is devoted to nonlinear finite element analysis of such 
apt The : aim 1 of the 2 analysis i is sto attain a apres insight i into the structural i 3 
The AXIPLANE ; program (13) is used for this purpose. In addition to the cracking 4 
of the concrete, this program considers the strain hardening and softening in 
the pre- and post-failure regions, respectively. The influence on the structural — 
= behavior of the uniaxial compressive strength, the ratio of tensile to compressive _ 
a strength, different failure criteria, and post- failure behaviors are ‘investigated. 
_ Moreover, because much dispute has been placed on the type of failure actually — 


occurring in the concrete, special attention is given to the structural behavior _ 


the failure mode. Finally, the predictions are compared with 


"Research Engr., , Risd National Lab., , Postbox 49, DK-4000, Roskilde, 
Note.—Discussion open n until September 1, 1981. Toe extend the 1e closing date one month, 
a written request must be filed with the Manager of Technical and Professional Publications, b 
ASCE. Manuscript was submitted for review for possible publication on April 3, 1980. | 
pi This paper is part of the Journal of the Structural Division, Proceedings of the American a 
Society of Civil ©ASCE, Vol. (107, No. ST4 1981. ISSN 0044- 
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‘ 
“applicable for axisymmetric and plane structures, is applied. Only axisymmetric 
triangular elements with linear displacement functions are utilized in the present 
In the numerical p performance, the total load is considered, even though it 
; is increased stepwise. However, for each load level, iterations are carried out 
until the constitutive equations for the concrete are in accordance with the 
- total loading in question. To consider the effect of stress redistribution, small 
load increments, around 2%-4% of the ultimate load, are employed. The — 
load is determined as that load in which a large number of iterations—i. e., 
25—is insufficient | to satisfy both the constitutive equations and ‘the 
‘equilibrium. 
_ The constitution model of the concrete is quite simple to ont 
is based on nonlinear elasticity, where the secant values of Young’s modulus 
and Poisson’s ratio are changed appropriately (11). This alternation is obtained . 
through use of a —— index th that relates: the actual Stress state to the 
Let Cas and o, o, denote the | principal stresses, and let o, > o, 2 
* which tensile stress is considered positive; then if only compressive stresses: pe 


present, the nonlinearity index, ‘B, is defined as fol. 


in which = ‘= the actual, most compressive principal stress; and oy the 


corresponding failure value, provided that the other principal stresses, c, and 
> are unchanged. Thus, 0 = B < 1, B = 1, and B > 1 correspond to ab 
states located inside, on, and the failure surface, ‘respectively. The 
"nonlinearity index, given by Eq. 1, is proportional to the stress for uniaxial 
compressive loading, and it can therefore be considered an effective — 
When tensile stresses are present, we transform the actual stress state ,, 
Zz o;), in which at least o, = a tensile stress, by superposing the hydrostatic _ 
pressure, —o,, and obtaining the new stress state (o/, 0/, 04) = (0, o, =o, 
a biaxial compressive stress We then define B as 


in which = the failure value of o}, provided that anda; are unchanged, 
i.e., the stress state (i, 3, is to satisfy the failure criterion. When tensile 


Stresses occur, the ‘nonlinearity index as defined by Eq. 2 is less than ~ a 


_ In Eqs. | and 2, a failure criterion is involved and in the AXIPLANE ae 
tw wo options exist: (1) The simple and well-known modified Coulomb criterion - 
proposed by Cowan (2); this criterion combines the Coulomb criterion—here 
with a friction angle equal to 37° and a maximum tensile stress cri criterion; and 
the more but also ‘more complicated criterion proposed 7 the 


on 00 onv parison of these two criteria” 
with e experimental data has been performed by ~ writer in Ref. 12. rE 
Letting the secant values of Young’s modulus and Poisson’s ratio depend 
on the nonlinearity index, the constitutive model predicts the strain hardening 
- before failure, and the failure itself and the strain softening in the post-failure 
. region. Smooth stress-strain curves are obtained, different post-failure behaviors — 


can easily be dealt with, and dilatation is simulated. Any failure criterion can 
_ be utilized, and the choice of an accurate criterion by itself assures accurate 
stress-strain curves. All stress states, including those where there are tensile | 
> _ inherent in the model. As described in Ref. 13, the applied implementation — 
of the model i in n the finite element program results i in a behavior 
to that: ofa 
compared to ‘that of strict nonlinear elasticity. ‘io. ‘de com 
_ The model is calibrated to a given concrete by six parameters only, which — 
are all derived from uniaxial tests. These parameters are: the initial mai 
i.e., the initial modulus Ey and the initial Poisson ratio, 


= parameter, D, that determines the degree of strain softening when the crushing — 


- Through the maximum tensile stress criterion, the modified Coulomb criteri 
contains in itself a cracking criterion given by o, => o,. o,. The failure criterio 


_of the writer (10) corresponds to a smooth surface in terms of one equation, 
_ and cracking is assumed to occur if: (1) The failure condition is violated; and 


(2) if o, = o,/2 holds. This cracking criterion was proposed in Ref. 11, and — 


its close agreement with experimental data is demonstrated in Ref. 12. The 
crack plane is assumed to be normal to the direction of the principal stress, 
q at the moment of cracking. Once a crack | is developed, it is assumed to 
remain open with a fixed direction. The standard smeared cracking approach 
is utilized. The shear stiffness along the crack plane is reduced to 4G, in which | 
Ge the shear modulus of the concrete; and y = the shear retention factor, 
0 =< y = 1. This factor is subject to much dispute, because in reality it is 


a a complicated function of crack-width, relative displacement tangential to the 4 


crack plane and the nature of the crack surface. For simplicity, however, we 
— consider the shear retention factor to be a fixed value, and y = 0.01 is utilized. 
The predictions obtained i in Ref. 13 suggest this value to be sufficiently accurate. 


The configuration of the specific pe out test considered here, the Lok-Test 
1. During application, 
a test bolt, consisting ~" a stem ond a circular steel ian, is mounted inside 
_ the form; see Fig. l(a). After curing the concrete, the form is stripped, and 
L - _ diameter than the stem is screwed into the disc and a cylindrical counter pressure 
is mounted; see Fig. 1(b). The rod is loaded by a pull-out force until failure, 


the stem is unscrewed. At the time of testing, a rod having a slightly smaller 


4 stresses can be dealt with, as shown in Ref. 11, but path independency is J 


4 the strain, €.(€. > 0), at uniaxial compeessive failure; and finally the post- 


+%. 


of 
| 


where. a a small piece of concrete « can n be punched out if: ‘ 
* the rod is applied. As shown in Figs. 1(b) and 2, this piece of concrete 
has the form of a frustrum of a cone. The meridians are almost straight lines — 


© “that connect the outer periphery of the disc with the inner periphery « of the 

_ Experimental data for the Lok-Test have shown a linear relation between — 

the force required to extract the embedded steel disc, and the uniaxial compressive 7 


- Strength of the concrete. A general status of the various oy rman investiga- 


‘mesh: consisting of 441 triangular elements. _ The elements that represent the 
steel disc appear in this figure. Perfect bond between the steel and the concrete 
_ is assumed. The pull-out force, as well as the boundary conditions at the location — 
_ of the cylindrical counter pressure, are also indicated. 


folure. 


FIG. 1 _—Application and a of Lok- Test (All Dimensions Arein millimeters) oo 


a To begin with, some i important aspects s of the structural behavior of the Lok-Test ‘ 7 
\ will be illustrated. Special attention will be given to the failure mode. After 
_ that, the influence of some concrete material data will be investigated in detail. — 
In the first place, the failure criterion of the writer (10) is utilized; later the 
effect of using the modified Coulomb criterion is evaluated. wre z 
To illustrate the structural behavior, we use concrete material data that can 
be considered as quite representative and realistic. For this purpose we approxi- 
mate the behavior of a specific concrete tested by Kupfer (7). The constitutive 
model of this relatively strong concrete is calibrated by the following parameters, 
all in accordance with experimental data: E, = 3.24 x 10° MPa; v, = 0.2; 
¢. = 31.8 MPa; o,/o, = 0. 10; «= 2.17%0; ond D = 0.2. Using these data, 
‘the normalized stress-strain curve is shown in Fig. 4. The values E = 2.05 
* 10° MPa, and v = 0.3 were employed for the steel disc. 
Let us first consider the predicted development of radial and circumferential 
cracks as the loading increases. This is shown in Fig. 5, where the loadings 
j are expressed in relation to the predicted failure load. The analysis ——. 
cracking in terms of cracked elements some 


3 
Bic 2 a the analuve tructure ac wall a the me ele art 
a 
3, 
4 


pl This type. of appears in Fig. 5(a), and is caused directly 
the pull-out force. At 18% loading, radial cracking initiates at the annulus near 
§ outer concrete surface. These cracks are caused by flexure similar to the 
- bending of plates. Such radial cracks appear in Fig. 5(b). With increased loading, 
a cracks shown there develop gradually. However, at 64% loading a considerable 
_ development of new circumferential cracks occur. As shown in Fig. S(c), wor 


concrete. Increased alien, | in particular, causes the radial cracks to develop; — 
the crack pattern just before — failure is shown in Fig. 5(@). pates ts 


FIG. 2. —Punched- Out Piece of Concrete (Cracks Are Made Visible somibesteaie Tracing) a 
While the predicted circumferential | cracks are supported by experimental 
‘evidence [compare Figs. 1(b) and 2], radial cracks have not been observed - 
experimentally prior to the present study. However, whereas the region theoreti- 
cally exposed to radial cracking is quite large, the corresponding crack widths 
are estimated to be quite small. If it is conservatively assumed that all tangential — 
deformation is concentrated in only one radial crack, then just before failure 


conclusion that no radial cracks are directly visible on the punched-out piece 

- of concrete. However, close inspection of such concrete specimens using r 

microscope reveals that clear radial cracking is indeed present. Such cracking | 
is visualized on the specimen in Fig. 2 using pencil tracing. 

w To further illustrate the structural behavior of the Lok-Test, the — 


| : are unloaded and crack “widths s reduced after failure, ‘thereby supporting the : 


a 
Gg the support. It is of interest to note that even though the concrete is severely ff 
7 ; cracked, its carrying capacity is far from being exhausted. Therefore, it is obvious 7 
| bs | 


distribution of 

_ the cracking is slightly more developed than indicated in Fig. 5(c). This stress . 
_ distribution is shown in Fig. 6, in which isostress curves are shown and the >) 


point. In accordance with the radial crack development, the distribution of the — 
tangential stresses shows large regions where tension exists. Only at the support, | 
and notably around the disc do compressive tangential stresses exist. The 
distribution of the maximum principal stress in the RZ plane also indicates 
_ large regions loaded in tension. Only in the vicinity of the disc, and notably 
/ at the support do small regions loaded in compression exist. The distribution 
of the minimum principal stress in the RZ plane is very interesting. Recalling 
2 


AND 


¢ 


that the uniaxial compressive strength of the concrete is 31.8 MPa, and noting 
that the loading is 70% of the predicted failure load, it appears that large stresses 


are present at the annulus near the disc. In fact, triaxial compression exists 


compression occasionally superposed by a small tensile stress appears in this 
region. Noting the stress directions, it is apparent that large compressive forces 
run from the disc in a rather narrow band towards the support, where triaxial 
as well as biaxial compression exist. This carrying mechanism is supported 
_ by the crack pattern of Figs. 5(c) and 5(d). It is of interest to note that both 


? “| __ directions of the principal stresses in the RZ plane are given at each nodal sy. 
| 
“4 
Moreover large compressive stresses are found at the outer — 
: _ of the steel disc, and comparison with the preceding figures shows that biaxial _ 
| 7 
circumferential cracks and the stress directions describe 


Ve 


PULL- OUT TEST 
a ‘slight curvature even ‘though they are almost straight. This small curvature 


is also observed in practice; compare to Fig. 2. works 
In conclusion, Fig. 6 shows that the stress distribution is very inhomogeneous. 
_ At increasing loading, a considerable redistribution of stresses can therefore © 
be expected. This suggests strain softening of the concrete in the post- failure 
region to be of importance for the failure load. Moreover, Fig. 6 shows that ' 


= compressive forces run from the disc in a rather narrow band towards» 


_ As demonstrated previously in Ref 13, the severity of the stress states is 
_ conveniently illustrated by means of the nonlinearity index, B, relating the actual — 

stress state to the failure surface; compare Eqs. | and 2. For composes 
loading, 0 = < |, = » and correspond to stress states located 
inside, on, an Fig. shows the 


ay 
3 


FIG. 4. —Normelized Strain tor @ (Stengths Are 


development « of the. increasing loading 
index, as a percentage. The distribution in Fig. 7(b) corresponds to the stress" 
distribution given in Fig. 6. Fig. 7 supports the preceding observations that 


loaded narrow band running from the outer periphery of the disc towards the 
support is very apparent. It should be recalled that when tensile stresses are 
_ present, the nonlinearity index is less than unity even at failure; compare to 

Eq. 2. At 64% loading, strain softening initiates below the steel disc, both — 
- adjacent to the annulus and at the outer periphery of the disc. At 79% loading, 
Strain softening develops from the outer periphery of the disc towards the suppo support. 

_ This development is pronounced at 88% and also at 100% loading; the latter 
= is the failure load, and the distribution given by Fig. 7(d) corresponds to the | 
last iteration before the calculations were terminated. At the failure _— 

| considerable strain softening occurs also at the disc adjacent to the annulus. 

can be as a decrease in n the index; ‘compare | Figs. 


00 -05 -10 -15 -20 -25 -30 -35 
| 
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the narrow region adjacent to the outer periphery of the disc, and running 
towards the support. This strain softening appears as a considerable drop of 
the nonlinearity index. This effect is very pronounced | when comparing Fig. 
-1(c) with Fig. 7(d), but a comparison of Fig. 7(b) with Fig. 7(c) already shows 
this tendency. It is important to realize that this gradual decrease c ofthe nonlinearity _ 
index due to strain softening in the post-failure region corresponds to the crushing — 
of the concrete. Thus, even though small tensile stresses may exist in addition 


and Hd). More important, however, is the strain softening occurring in 


FIG. 5.—Crack Development with Increasing Loadings (Loading Is Expressed in 
‘ Relation to Predicted Failure Load): fa) Loading = 15%; (b) Loading = (25%; (ec) 
Loading = 64%; (d) Loading = 98% 
to the primary biaxial compressive stress states, the failure is caused by the 


7 crushing of the concrete and not by cracking. This mechanism is in accordance ~ 
— the experimentally observed ductile failure mode. . Therefore, the force — 


may have some direct influence, ‘as will be examined later on. 
_ Let us now compare the predicted failure load with experimental data. Based 
on the results of different test series, including a total of 1,100 Lok-Tests, 


Kischegnesd-Hensen a Hansen and Bickley ( (6) ‘suggest the the + following | linear relation t between ” 


| 


PULL- -OUT TEST vest 
oie out F, ‘uniaxial ‘compressive cylinder swength; 
0.8 o., in which F and o. are measured in kilonewtons and segepascaa, 
respectively. This relation is shown in Fig. 8, and is based on concrete mixes, 
in which o, ranges from 6 MPa-53 MPa. The failure load resulting from the — 
™ - yieeate calculation, in which o, = 31.8 MPa, is also indicated. The analysis 
underestimates the experimental failure load by only 1%. 
To investigate the dependence of the oe value, a calculation was performed 
with data from another, weaker concrete. To wae use of realistic concrete 
Oe Wns heving he | 
Eves fat 


max. principal stress in RZ - plane 


tension 


cipal st RZ- pl 
FIG. 6.—lsostress Curves and Directions of Three Pri Principal Stresses at 70% Loading — 

(Quantities Are in megapascals) ol 

- data, test results of Kupfer (7) we were utilized again. In the constitutive model, 
following parameters were applied: E, = 2.89 x 10* MPa; v, 0.19; 


7 = 18.7 MPa; o,/o,. = 0.10; €. = 1.87 %o; on D = 0.6. The close agreement 
of the resulting predictions with the experimental data of Kupfer (7) has previously 
; been demonstrated in Ref. 11. In general, the weenee the concrete, the more 
ductile i is its post- failure behavior; compare to, €.8., et This 


a 
q 


ars? using these concrete parameters underestimates the ‘actual failure load iy _ 
only 3%, and is plotted in Fig. 8. Therefore, the calculations are in oe ena 

_ with the experimental evidence, showing that within the considered variation 

of the o. values, a linear relation exists between pull-c - fores and and compressive 

~ Itis remarkable that the prolongation of the: experimental line in Fig. 8 intersects 
the ordinate axis at some distance the origin. two of 


q 


FIG. 7.—Development of Contour Lines for | asa 
- (Loadings Expressed as a Percentage of Predicted Failure Load): (a) Loading = 20%; 


(b) Loading = 70%; (c) Loading = 94%; (d) Loading = 100% | 


on 

concrete are dependent on compressive strength, namely the ductility 
and the ratio of tensile to compressive strength. As has already been touched 
upon, the weaker the concrete, the more ductile the post-failure behavior. To 

_ investigate the influence of minor variations in the post-failure behavior of the 

concrete, a calculation was s performed u using again the concrete having a a strength i 

of 18.7 MPa, and possessing unchanged properties except for a lesser ductility. 

Thesefore, the value ‘D = 0 was used instead of the more realistic one, D 


| 

30 

| 

_ .. 5%, ‘as shown shown in Fig. 8. That the fi failure load depends on the particular fim 


softening behavior of the concrete is is indeed | to be | expected, ¢ ee the 
a general, the weaker the concrete, the larger is the ratio of tensile to 
compressive strength; compare to, e. g., Wastiels (16). Let us investigate this 
Si effect using the concrete having o, - 18.7 MPa, and D = 0.6 again, but putting 
H now o,/a, = 0.12, instead of o, / he 0.10. This increases the predicted failure 
_ load by 11%, as shown in Fig. 8. In reality, Kupfer (7) determined the o Jo. 
ratio to be 0.105 for the concrete considered, and if interpolation is performed — : 
‘between the two calculations having the o, [Oe ratio equal to 0.10 and 0. * 
7 respectively, then the resulting failure load is 0.7% below the actual value. 
Even though the tensile strength of the concrete certainly has an influence 
on the failure load of a Lok-Test, it is of importance to realize that this influence 
is an indirect one. Only very little of the pull-o out force is carried directly by 


Hansen and (1978) — 


D 


2/0. 10 


«id mod. d.Coulomb 
theoretical, sen ‘and Brestrup (1976) 
= 


8. Data Compared with Theoretion! Failure Values 


tension in tiie concrete, but the regions where failure takes place are are primarily 
s biaxial compression, occasionally superposed by small tensile stress. The rE 
a failure is caused by crushing, but the existence of even a small tensile stress" 


considerably decreases the failure strength, | 
‘The analysis has demonstrated that the reason that the relation 


a result. of the increasing ductility and the increasing ratio of tensile to” 


compressive strength, the weaker the concrete. 
Let us now investigate the influence of different failure criteria. For this ? 
_ purpose we return to the concrete having o. = 31.8 MPa. All constitutive 
_ parameters are unchanged, but now the modified Coulomb criterion is applied. — 


_ Compared to the previous analysis, this reduces the predicted failure load by 


| 
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in biaxial the modified Coulomb criterion is | 

to underestimate the failure stresses for such stress states by 25%- -30%. a 

It is of interest to observe that the decrease of failure load, when | using 
the modified Coulomb criterion, is an accordance with the finding that the 
_ Lok- Test depends directly on the compressive strength of the concrete and 
‘not on its tensile strength. As demonstrated in Ref. 12, the modified Coulomb 
. in general, underestimates the failure stresses, when concrete is loaded — 

ie compression, except when extremely large triaxial compressive Stress states 
exist. Moreover, as shown there, this criterion, in general, overestimates the ‘ 
failure stresses when tensile stresses are present. Therefore, if the failure in . . 

a Lok-Test was caused by tensile cracking, then use of the modified Coulomb | yr 
criterion would result in an increased failure load. However, in accordance 
with the preceding analysis, use of the modified Coulomb criterion decreases 


Jensen and Brazstrup (4) have determined the failure load 
a Lok-Test using rigid-ideal plasticity theory. _ They also used the modified 


ii Coulomb criterion, and their result is shown in Fig. 8. It appears that in 


_ agreement is obtained even though proportionality and not just linearity cell 


a: the pull-out force and the compressive strength was obtained. However, the tsa 


failure load determined by Jensen and Brazstrup (4), when o. =31 8 MPa, 
is considerably larger than the one determined here when using the modified 
= criterion also. This is particularly conspicuous, because Jensen and 
Beyer y (4), in their analysis, are forced to use a friction angle equal oS 
e angle « shown in Fig. 1(b). This results in a friction angle, @ = 31°. The 
_ present finite element analysis is based on the value @ = 37° which, as previously — 
examined, results in some underestimate of the actual failure stresses. Use 
of the value = 31° would indeed infer a considerable underestimate of actual 
- failure stresses. However, Jensen and Brazstrup (4) in reality compensate for 
: this, as their analysis is based on rigid-ideal plasticity with no softening effects 
at all. Consequently, they assume failure all along the plane running from the - 
Bek periphery of the disc towards the inner periphery of the support. Previous 


examinations related to Fig. 7 have . refuted such an assumption. However, in 
. accordance with findings in Ref. 13, this underlines the extreme importance 


of including a suitable strain softening behavior in ramen modeling of 


In the present takes of Ref. 12, the structural behavior 

of a specific pull-out test, ie., the Lok-Test, has been investigated in detail. ue 
: Severe cracking occurs, and the stress distribution is very inhomogeneous. It fs 
is has been shown that large compressive forces run from the disc in a rather 
= narrow band towards the support, and this constitutes the load-carrying ie 7 

nism. Moreover, the failure in a Lok-Test is caused by the crushing the concrete er 
ps and not by cracking. Therefore, the force required to extract the embedded 
7 steel disc is directly dependent on the compressive strength of the concrete 

in question. However, as the stress states, in which failure takes place, are 
Primarily | biaxial compressive occasionally superposed by tensile ox. 
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‘the tensile “strength of the concrete has some influence. ‘The effect 
of strain softening in the post-failure region is important. In general, — 
concrete compared to strong concrete has a relatively larger tensile strength 
and a higher ductility. This explains why the relation between the failure pull- -out - 
_ force and the compressive strength is linear and not proportional. cits 
a With | respect to finite element modeling, it has been demonstrated that not — 
— only is use of an accurate failure criterion important, but modeling of the strain _ 


_ softening in the post-failure region turns out to be a mandatory prerequisite a 
for realistic structural predictions. 
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By Saleeb' and W. F. Chen,? M. ASCE. 


aa Large displacement analysis of elastic- plastic pipe bending problems is fre- 
quently encountered in many offshore structural engineering applications. One 
ss such application is the bend around of pipelines during a pull-in operation. — 
In this case, a monotonically-increasing pulling force is applied at the tip of 


_ the pipeline in the sea, and is directed toward a fixed pull-i in base | point of 
offshore platform (located by the so-called “‘overlength and offset”’ distances, 
_as shown in Fig. 1). The horizontal component of the pulling force (acting — 
parallel to the pei pipe position) will produce a large secondary bending @ 
moment (commonly called the P-A moment) in addition to the primary bending — ‘ 
moment due to the vertical component of the pulling force (which : acts perpen- | 
dicular to the original pipe position). Thus, the additional P-A moment must 
_ be included in the analysis of the pipe bending problem. & ot 

_ It is obvious from the nature of the problem that the required final ha 


- _ deflection is very large (see Fig. 1). Therefore, the more exact expression wil 


“ the curvature, ®, and deflection, ,, may have to be employed in the analysis. 
such a problem. requires the use of the exact but 


instead of 
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- 1.—Plan View of w of Pipeline Being. Pulled into Fixed Base and Definition of Offset 


in which R = the radius of curvature; we dw/dx; 
= the distance measured along the center line of the 
_ To resort to formal mathematics in solving such a problem requires 2 a rather 
~ complicated analysis which always | leads to a nonlinear differential equation 
coupled with the displacement-dependent lateral frictional forces due to | pipe-soil — 
interaction. This nonlinear differential equation is often intractable even by 
_ the well-known numerical technique of the finite difference method. A more 
- elegant approach would be to make physical idealizations of the pipe, coupled 


With | a ‘numerical procedure, in order obtain solutions for practical 1 ‘This 


APPROACH AI NSIDERATIONS 
‘The problem 1 of pipe bending analysis will be iaiinaaiiial’ from the por 
of beam-column theory with large displacement analysis, using the finite segment | 
method and incremental tangent stiffness procedure (see chapter 11 of Ref. 
6. More rigorous approaches, such as shell type of analysis is not considered 
a here. It is believed that the pipe bending problem is essentially one-dimensional, 
_and therefore can be best approached from this type of analysis. a Oo 
_ The present work covers the in-plane behavior of the pipe under static loadings 
only. . Pipe dynamics during pull-in will not be considered. The stress-strain | 
relationship of the steel is assumed to be of the elastic-perfectly plastic type. — 
In addition, torsional effects on the pipe tubular section are ignored because ~ 
_ of the large torsional rigidity of the closed thin-walled tubular section. 
AS for the pipe-soil interaction model, only a simplified approach “will be 
Be a: here. In this approach, the pipe-soil interaction is simplified as a series" 
of lumped concentrated friction loads at the nodal points along pipe length; 
- a of these point friction loads follows the Coulomb friction relation with ss 
single value of friction coefficient for both static and Sliding friction. 
_ The beam-column analysis of the pipe bending problem, as presented here, q 
| . can be divided into three steps: (1) The moment- curvature behavior of a short 
Segment, (2) the pipe-soil interaction modeling; and (3) the finite segment — a 
formulation, i.e., the segment stiffness matrices. These will be described briefly 
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a Moment- Curvature (M- ®) Relations. —The basic quantity pied in any a 
— beam-column analysis is the value EJ, which here can be considered as the 
slope of the relation curve between moment, M, and curvature, ®. In the elastic | 
sone this quantity is constant and therefore presents no difficulties. However, — 

_ in the plastic range, this quantity is not constant since the moment-curvature i 
oy is nonlinear, and the instantaneous bending rigidity, EI, should be — 
used. Even if the stress-strain relationship of the material is assumed to be 
an elastic- -perfectly plastic type [see Fig. 2(a)], the generalized stress-strain ts 
_ M-® relationship of the segment is seen to show an elastic strain hardening — 
characteristic, as shown in Fig. 2(b). This is due to the fact that the section 
_ can be partially elastic and partially plastic. 


idealized Strain Relationship (Elastic-Perfectly Plastic Typical Mo- 


& In the case of pipeline of fabricated ‘tubular cr cross section, atus are many &. 
factors which influence the M-® relations. These may be summarized in the i 


following (see a 6 of Ref. 6; see also Refs. 7 and 9): (1) The intensity 
of axial load, ee in which 1 P, = the yield axial load | of the cross “section; 


the process - fabrication; (3) the initial out- of-roundness; and (4) the diameter 


wT. 


__ Inorder to obtain the moment-curvature relations which include the interactions 


between these factors, a a a computer subroutine has been developed and use 


aan 1981 T4 
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inRef.9 (using the tangent stiffness method). However, the use of this etieutinn 


_ A brief description of the method outlined by the ‘second writer (3) for obtaining 
approximate expressions for the moment-curvature-thrust (M-®-P) relations, — 
and their specific applications to fabricated tubular cross section, are given — 

_ The general M-®-P relation as shown in Fig. 3 can be divided into three 
regimes: (1) Elastic; (2) primary secondary plastic. Using the 
_ nondimensionalized quantities, p = P/P, = M/M,, and 6 = 9/®, ~ 


which P,, M,, and o, = the yield load, vield moment, and yield curvature 


consiony 
\secondory 


3. —Moment- Curvature Thrust Relationship for Common Structural Section 
« Bis 
‘% of the cross s section, respectively, the not nonlinear moment-curvature relation for 


elastic-plastic sections may be represented b 
4 ee = =o pike 


3 


m= 
The parameters a, b, c, f, m,., b,, and , are functions of thrust, p, which cho 
a are determined by test results or numerical computation for a ‘particular shape pe 
of cross section. Once the values of m,, m,, m,., &,, and , are known, . 
an = parameters a, b, c, and f can easily be determined by solving the two - 
sets of two simultaneous equations which arise if the particular values of >, 


and b, are inserted in the appropriate equations, i.e., , Eq. 3, and the moments . 
4 —_ to the appropriate moments, im, and m,. The results are given by © 


1 
‘usually results in the most time-consuming part Of the Caiculations. computa- of 
q tional time will be reduced significantly if closed-form moment-curvature expres- ; 
a 
4 
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Approximate m-d-p Expressions for Tubular Section.—Exact analytical ex- 
pressions of m,, m,, m,., ,, and ¢, for a tubular section subjected to axial 


load, P, and bending moment, M, have been developed, and used for plotting * 
: 4 the interaction curves of ™ 4, These expressions are very complicated and 


FIG. 4.—Interaction Curves for Tubular 


> + 1.206 
for 50.4; 1.528 (1 2.05 (1 - 


= 1.82(1 82(1-p), for 0.65.5 
As an example of m--p relationships obtained using Eqs. 
5b, and Sc, together with Eqs. 3 and 4, the approximate m-¢-p curves for 
_ the tubular cross section for four cases (p = 0, 0.4, 0.6, 0.8) are plotted in— 
‘Fig. 5, in which | actual m--p relations, as reported in Ref. 9, are also included 


| 
be used to approximate the exact interaction Curves (12). 
| 
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Actual (Chen/Toma 1979) 


. It is seen that the curves are” almost identical for both the 
Approximate m-¢-p Expressions Including Effect of Residual Stresses.—Based 
on the results obtained in Ref. 9, the approximate m-p-d relations, including 
the residual stresses effect, can be derived. ‘The details of the longitudinal and 


= 0.893 — 1.305 p + 0.796 p’, = 0.893 1.305 p +0.796p, 
“for 0=p<06; m, = 1.50 2.333 p + 0.833 9°, >, 


0.833 for 0.6= p= 


= 0.893 + 0.546 p 1.946 p?, 

for O<p=06; m,=1 1,375 - 150p +0.125p?,6,= 1.040.875 

& 

1.273(1 — 1.18p p< 0.65; ; m,.=1 


| 
per ,, and can be closely approximated 
= 
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AIG. 6 Thrust Relations (with 
e curves s curves are re shown i in Fig. 6, together with the ecteal relations| 
- = (9). Again, the approximate m-¢-p relations agree well with the actual relatio ion- 
om Finally, similar approximate | m-o-p expressions, including the effects of residual i 
_ stresses, out-of-roundness, and hydrostatic pressures can be developed. This 
i is reported elsewhere (9,12), Be 
_ Pipe-Soil Interaction Modeling. —As mentioned earlier, the pipe-soil interaction 
will be simplified as point friction nodal forces, each of which follows Coulomb 
friction relation. Although the coefficient of friction will have different static 
dynamic values, , a single value for both static and sliding friction will be 


The preceding is, of course, an extremely : simplified approach for modeling _ 


_ the pipe-soil interaction behavior; nevertheless the general rigid-plastic type of — 
load deflection behavior is valid and acceptable. In order to follow the actual — 


load deflection behavior, it is necessary to conduct extensive experimental s studies — 
which will give the complete load deflection response from zero drag force - 
upwards to sliding force, taking account of various directions of the pulling 7 
force. Thus, if this simplified drag force displacement relation is not acceptable, 
then a nonlinear or a piecewise linear relation must be developed using curve a 
fitting techniques with experimental data, coupled with some simplified rigid- _ 
. plastic type of analysis, such as limit analysis of soil plasticity (4). This is 7 
not considered inthe present work, ab 
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approech adopted here uses a set of constraints to impose approximately 
the rigid no sliding condition, or the constant force sliding conditions. In order 
to approximately model these conditions, a spring element with a relatively 
high stiffness coefficient is used. In this model, the load displacement relation 
of the spring is of the linearly | elastic- perfectly plastic type, as shown by Fig. 
7. As long as the friction force does not exceed the friction limit, F,, the 
4 spring has a constant stiffness coefficient, k,. Thus, the rigid no sliding (stick) 
_ condition is not imposed exactly, but approximated with an acceptable error. 
When the friction force exceeds its limit, the spring stiffness becomes zero 
and slip occurs, i.e., , large displacements occur without further increase in the 
“friction force. A refined alternative approach for implementing the friction <_< oe 
oF slip conditions can be found in Ref. 10. 
a Stiffness Coefficient of Spring Model.—The choice of the spring ‘stiffness 
coefficient is extremely important in the pipe bending analysis, since it has 
been found that the convergence rate of the finite et incremental solution © 


Load (F) 


Bostic (approximate stick condition) 


7.—Simplfied Load Relation for Elastic- Perfectly Plastic Spring 


Used for Friction Modeling _ 


“presented here is very sensitive to this choice. in fact, using too high ssl 


coefficients causes the spring elements to dominate the solution excessively 
and produces a very large jump at slip. On the other hand, the use of too 
low stiffness causes an appreciable error in modeling the rigid no ) sliding (stick) — y 


condition. Therefore, a compromise is essential between accurate modeling of ] 
the stick condition and avoiding spurious effects of a too dominant stiffness _ 
and large jump at slip. Thus, the argument that has been used is to determine _ nf 
an acceptable tolerance in the stick condition approximation. - Referring to to Fig. — 
7, the the stiffness coefficient, + in the linear elastic range is is given t by lie dale 


in which p = the friction coefficient; N = total normal force; and Din = 
the ‘‘assumed”’ displacement limit. This displacement I limit mnteugueeeaeee the assumed 
error in modeling the stick condition. — 
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More Refined Friction Spring M Model ~The in of the friction 
e force ; during pull-in of the pipeline is not considered in the previous model. ., 


To account for this change, a more refined spring model may be used. In — \ 


_ this case, a spring model having two stiffness components is used, as shown 
_ in Fig. 8. The resultant frictional force is assumed to bei in a direction perpendicular 
to the pipe segment at the nodal points. Therefore, , the stiffness components, — 
ki and ki, in Fig. 8 are: k, = k’ cos ; ki, = k’ sin y, in which k’ = the 
‘resultant stiffness coefficient which has the value ky (see Eq. 7) in the linear 
_ elastic range; and W = the notation angle of the pipe segment as oe in 


Finite Segment Formulation.— = the finite segment method, the actual tubular - 
pipeline i is physically replaced | by an assembly of finite segments, and the problem — 
is formulated and solved approximately i in terms of the behavior of these segments , 
_ without recourse to complex differential equations. In this sense, the finite 
_ segment approach may be considered as a simplification of the more general _ 


of the finiteelement method (14), 


In the analysis of pipe bending problems by the finite segment method, — 
and rotations are into account at each nodal cross" section. 


friction 
direction 


‘Stiffness 


Component 


In computation of the segment matrices, the direct stiffness 
is utilized. Geometric and material nonlinearities are handled by application 


of a modified tangent stiffness apr approach. The details of the method and its” 
~ application to solve biaxially-loaded beam-columns of — — shapes in 


and o one rotation a are taken into a account at each nodal cross section (see Fig. > 


+ 9). The fellowing assumptions are made in developing the segment stiffness 


_ behavior; (2) section properties for the tangent stiffness of a segment are invariants 
4 within the segment, and a segment tangent stiffness is calculated from the slope 
of the M-® relationship; (3) twisting effects are ignored; and (4) material 

Stress-strain relation is linear elastic-perfectly plastic. 
Cm _ Segment Stiffness Matrices.—Consider an undeformed segment, AB, of length 


Las shown in Fig. 9, with the local coordinate axes as the 


oblems, the pipe will be al 
4 


ina 
principal axes of the elastic portion of the cross section _ At each nodal cross 


section, there are three nodal forces and three nodal displacements, a 
m,, and u,, u,, respectively . The segment stiffness | matrix: [K(6 
6)] relates the six nodal forces, {f(6)}, to the 
in which {f} = the ‘nodal fence vector of components Rev 


and an and {u} = = the nodal displacement vector of of components u,,, 


ft 


sey segment “stiffness ‘matrix [K(6 x is ‘computed, considering axial 


FIG. 9. —F Displacements of (Positive Directions for Forces and 


deformation and bending deformation separately. Further, the section properties 


are assumed to be constant within the segment. Thus, [K(6 x 6)] is derived — 
~ by solving, separately, the axial and ae problems, using the following 


i 
| 
| 
4 
which the primes denote differentiation witn respec equation 
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FIG. 10. of (a) Axial Deformations; 
. 10 governs the problem of an axially-stressed segment, Fig. 10(a), while x ; 
_ the second equation governs the flexure of a beam-column Segment, ‘Fig. 10(b). 
2. A set of six boundary conditions (three for each terminal cross section) should . 
be attached to Eqs. 10 in order to define the boundary value problem. The 
_ final results obtained are summarized 
Different forms of ‘submatrices. [Kal [K,.]. and are are given 
in the following for various axial forces (compressive, tensile, and small): 


El, k(s —kLe) 


oh: 


ee 
| 


4 
& 
‘ 
4 
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—- 2. For Axial Tensile Force: (P = fea * =f, 


APRIL 
antisymmetric 


El -c) 


in which k= P/ EI, cos kL: s = sin kL; end of inertia > 
of the cross section about the €-axis. MOA, 
ce) 


(0) 


> 


| 


antisymmetric 


“* 


A = cosh kL: and s’ = kL. 


For Small Axial Force (Tensile or The factor 
— 2c — kLs) or (2 — 2c’ + kLs’) in the previous stiffness matrices approaches 
“zero as the axial force becomes very small. In order to avoid numerical errors - 
in such cases, an explicit form for the stiffness matrix is . oe ‘Using — 


Taylor ser 


for small values 


ta | 
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kL), en stiffness matrix ix reduces to familiar form a beam segment, 


laboa 

<_< <= EA 

6EI, 


_ The basic feature of a problem having geometric seatndaitiny' is that equilibrium 


ee equations must be written with respect to the deformed configuration, | which © 
is not known in advance. Therefore, the element-to-structure transformations 
_ are based on the deformed geometry of the element. Referring to Fig. 11, the - 
 Jocal and global coordinate axes are &, n, ¢, and X, Y, Z, respectively. After 
the global total displacements, U Un, and U,,, have taken place, 


the undeformed A 0 Bo, assumes | the position AB (see 


sd 

af) suk fur 3 


11 —Displecements of Segment. (in- -Plane Behavior) tom wit 
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in which vf is positive if it is measured clockwise at the left end, B, of the 


"Segment. ‘The global nodal forces, {F(6)}, and displacements, {U(6)}, are 
to local vectors, {f} and {u}, through the relations: 


(U) =[R}"(u); and 

in which [R(6 x ball “ = the — of the rotation matrix, [RO x 6), 


hich is gi | 
_ whic is given by: a’ 


=| [RG x 100 x: x 


[RO x 6)] 
(03 x3) [RG x3)] 


which siny cosy 


[0(3 x 3)] is a3 x 3 null 


oat 
value of the segment bending stiffness, EJ, of the cross section can 
( dl be easily determined from the slope of the e nondimensionalized moment-curvature — 
7 relations (by differentiating these relations which respect to ‘the nondimensiona- 
lized curvature, 6). However, the segment axial stiffness, EA, must be determined 
a by considering the elastic part of the cross section only (for linear elastic-perfectly 
plastic material). After ag -_r of the cross section, the instantaneous 


in which 7 and Io = the | area and the moment of inertia ‘about | the oe 
of the total cross section. The variables, B, and B,, may be considered as 


reduction factors for axial and bending stiffnesses, respectively, due to the 
_plastification of the cross section. For thin-walled tubular section under pure 


3 bending, it can be shown that for B, and B, are 


1 foe for $= 


Variables 8 , Ba + | (1/7) sin (2 cos" 1/o); and B, = 


(1/%) — (1/7) sin (2 cos~' 1/6), for to 


= tan” ' ———— 
— 
| 
&§ 

S(N x N)] =SIR] IKIIR 
7 
q 
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ANALYSIS OF PIPES 


IcnemenTaL Mooirien Tancent StirNess mess tas bes 
7 ‘The problem to be solved now may be ony as follows: Given the state, : 
Ay a at which a ‘structure is subjected to forces {F, }, and has displacements. : 
4 {U, }, and tangent stiffness [S,], as shown in Fig. 12, it is required to find 
; a the displacements {U,}, when forces are increased to {F,}. Because of the 
7 nonlinearity of the system, the correct values of {U, } cannot be obtained directly. 
Herein a modified tangent stiffness approach is used (Ref. 6, chapter 11). The 
stiffness of the system during any loading increment is assumed to ) decrease — 


monotonically. Therefore, the secant stiffness, between een the two s states 
ar A and B, is bounded by |S,| = 5 = 


" {Fs }, is summarized i in the following steps (see Fig. 12). ner q diet) on 


iff 


| 
a for sma 
j 


Using t the initial unbalanced force vector, {AF, {F, _(F, 


: the current tangent stiffness at A, [S,], «], solve for the incremental displacement, _ 

“Calculate the new displacements, {U,} = {U,} + {AU, 4 

_ 3. From the updated displacements, {U,}, the updated tangent stiffness, [S,], 


wee 


620, 


” ie _ Repeat the same procedure to obtain the new state represented by point 

in Fig. 12, i.e., {AU,} = [S,] “'{AF,}, (F,} = [a S, + (1 
.. {F,}, and {AF,} = {F, }- {F,}. Thus, after n cycles of iteration 

ag cycle, n, check the error, (AF, }. Stopif the unbalanced 

a 


ecome less than a prescribed tolerable limit. egal tae 5 


oral 
* Based upon the tibiinatoad formulation, a computer program has been 
developed at Purdue University to provide numerical results. A number of —_ 


_ whenever possible, with the available theoretically obtained solutions. Finally, 
as a practical case, the developed computer program is used to carry out the — 
‘ elastic analysis for a large displacement pipe bending problem. Some typical 

Check Problems.—The results | for comme of of the check problems a are 

‘summarized in the following: 


2 Large deflection analysis of an elastic cantilever beam: The exact solution ny 
of this problem i is obtained by Bisshopp and Drucker (2), where t both the exact 
a curvature expression and the effect of shortening the moment arm as s the loaded 


end of the beam are The same is solved using 


a? 
| q 
(check) problems are solved to provide the necessary confirmation of the validity 
. _ of the proposed numerical procedure. The results obtained are compared, 7 
Beam: Effect of Load Steps 


—a=05(4 (4 segments) 


Exact solution ( Atsuta, ome 


Case(3) Load P at L/6 from en 


FIG. 15. —Beam or on Elastic Supports with One Lateral side Load at eat 7 


oh the assumption Pre a smail displacement derivative during each load mernen, 


hich is the basis of deriving the governing differential equations (segment 


ANALYSIS OF PIPES 

‘ _ the developed computer program. Typical results are shown in Fig. 13 | 
vertical and horizontal displacements, 6 and A, respectively, of the free end. =f 

| As can be seen, the finite segment solution agrees well with the theoretical a 

results (2). The ion 

the ses 
— 
FIG. 14.—Pin-Ended Column with Initial Deflection (Effect of Averaging Factor a) 

Load Pat midspon L/2 — 

— | 


tions (out-of-straightness): Typical results are shown in Fig. 14, where the exact 
- solution of Ref. 1 is also plotted. A good agreement is observed. _ oc SA 


“Load at 


‘Load ot L/3 


Exact ( Hetenyi, 946) 


FIG. 16. —Deflection enone for Beam ‘on Elastic Supports with Different. Load 


M, 


segments 
(8 segments, a= 0.5) a 
— 00050010 
rotation 6, (Radians) 
FIG. 17.—End Moment-End Rotation Relation for en Column of Rectangular Cross 
Section Subjected to End Moments 
_ Finite beam on elastic supports: For different locations of the load, P, 
the beam (see Fi the calculated « curves are given in Fig. 


End moments 


| 
Be 
Cross section: Rectangular D = 0.06L mn? = 
| 
ta 


= 


16. Finite segment results are seen to be identical to those of the theoretical — 


on 4. Elastic- -plastic analysis of beam-column with rectangular cross sections: ns: 
od 


- The finite segments solutions and the results obtained using Newmark’s meth 
(5) are plotted in Fig. 17. Again, a good agreement is observed. 
In short, the numerical results pevcented: demonstrate the validity of the 
computer model i in solving dif types of having mend 


the s sea floor. The s ‘descriptive data are | given in Fig. (18. The 


is modeled using the mesh shown in pan 19. The tubing-floor interaction is 


- which are modeled by the springs shown j in Fig. 19 (tubing weight = 19.5 lb/ft 
_ or 284.5 N/m). The stiffness coefficient of the spring models is determined 
a on the basis of an assumed value of fiw | equal to 3 in. (76.2 2 mm). ‘a 

ty Overlength of 


tine 


Inside diameter D;= 4.276 in. ‘ ae 
stress oy = 90 ksi 


fe. 18.—Top View of Tubing Pulled In on Floor: Dimensions and Properties a 


mm; 1 ksi = 6.89 MPa) ot 24 


an Number of nodes = 27 hares 


| 
| 
.—l 
g 
q 


APRIL 


‘spring , models are . considered: ‘the ‘first model with one spring stiffness ii in the 


vertical direction, and the second model with two spring stiffnesses in the vertical a 


and horizontal directions; each has the same friction force limit (extreme case 


a 


aw 


— 30 90120 150180 


ME 


Two spring stiffness 


Tip rotation angle @ 


7 trials were made before the final choice of this value. = = = = | L 
n for Tubing (Overlength | 
Offset = 190 ft. 
pit 
q 
19004 
— 
FIG. 21.—Pull-In Force-Tip Rotatio 
1 ft = 0.305 m;1 Ib = 4.45N) 
| 


sm 


Fina rotation @= = 92° 


II: Configuration at intermediate pull-in force ( pull-in= 3000 pounds) 


Final configuration at end of pull-in (pull-in= 4200 pounds) 
FIG. 22.—Tubing Configuration During P Pull-in (Offect = = 190 ft; 0 Overtength = 
is Typical force- -displacement responses of tubing obtained from the finite segment 
‘solution are shown in Figs. 20 and 21 . Tubing configuration at final pull-i in 
’ and at an . intermediate value of pull-i in force is given in Fig. 22. > Shel dale bey 
_ solutions obtained for all the cases investigated indicate that friction F 
modeling has noticeable effects on the results. Numerical difficulties are — 
< tered in some cases and a large number of loading steps are used. Therefore, 
q great deal of work and further studies are needed in this matter to achieve 


a reasonable specific conclusion concerning. this effect i in | ‘the he analysis, of actual 


Large displacement analysis of pipe bending problems i is frequently encountered 
in offshore structural applications. Any r rigorous analysis which "attempts to 
solve this type of problems is destined to be complex, since both nonlinearity — m4 
~ due to material plasticity and nonlinearity due to geometrical change must be | 
included in such a formulation. Herein, the numerical procedure based on finite | 
segment method is presented. The pipe is treated as an assemblage of ‘‘cans’”’ 
‘segments’’ and the moment- -curvature Telation of the tube cross 
is approximated by simple mathematical expressions. ‘The segment stiffness 

relationships are computed using the direct stiffness method. Material and — 
geometry nonlinearities are overcome by a modified tangent stiffness approach. 
_ Four numerical check problems are then presented and compared with “‘exact a 
_ solutions. A good agreement is observed for all cases studied. Finally, the problem 
of bend-around of pipelines during a pull-in operation is studied in detail. ol 


_ method is found to be adequate and efficient for computer solutions. =~ 
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INTERMEDIATE STIFFENE RS FOR THIN-WALLED is 


By Thomas F Desmond,’ "A. M. ASCE, Pekoz,’ M. age ef 
George Winter,” Hon. M. ASCE 


This paper is of a two paper sequence describing the of 
an investigation of the structural behavior of longitudinally stiffened compression 
elements of thin-walled members. The first paper (5) dealt with edge- “stiffened 
- elements, while the present paper deals with intermediately stiffened elements. 79 
In both cases, the ‘stiffener is longitudinal, i.e., parallel to” the direction = 
- the applied in- -plane c compressive stresses. A detailed description of the complete 
investigation may be foundin 
A typical component element of a thin-walled member is the stiffened element 
4 _ shown in Fig. l(a). A stiffened plate element by definition is adequately supported © 


on each longitudinal edge. Local plate buckling is a usual design consideration 


for thin- -walled members. One way | to increase the local buckling load and to 
enhance structural efficiency i is by providing a secondary or intermediate stiffener Pod 
- between the main longitudinal stiffeners or webs. As shown in Fig. 1(b), the © . 
_ intermediate stiffener breaks up the wave- “like pattern of the stiffened pa tear 
_ of Fig. 1(a@) such that the two plate strips, one on either side of the intermediate 
_Stiffener, behave as two stiffened elements. 


a 


intermediately stiffened elements analytically and experimentally. Based suivante » 
on the experimental results, stiffener requirements were developed. Methods _ 
were also formulated for predicting ultimate strengths of elements that are 


supported by such stiffeners. = 


_|Formerly, Research Asst., Cornell Univ.,Ithaca,N.Y. 
? Assoc. Prof. of Structural Engrg., Cornell Univ., School of Civ. and Environmental 


Engrg., Hollister Hall, Ithaca, N.Y. 14853. 
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7. Note.—Discussion open until September 1, 1981. To extend the closing date one month, 
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ASCE. Manuscript was submitted for review for possible publication on April 30, i 
This paper is part of the Journal of the Structural Division, Proceedings of the American Pi 
Society of Civil Engineers, _@ASCE, Vol. 107, No. ST4, April, 1981. ISSN 0044- 
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Two distinct buckling modes characterize the buckling behavior of interme- 
diately stiffened elements. One is the stiffener buckling mode, where the instability — 
‘ of the assembly is initiated by buckling of the flange-stiffener assembly direction “i 
perpendicular to the plane of the flange. Stiffener buckling will ay 
induce local buckling of the component plate elements. The second mode is 
— the local plate buckling mode. Both buckling modes modes are shown schematically a 
_ buckling mode. For the stiffener buckling mode the buckling equation is given 
by Barbré (2), and for the local plate buckling mode the equation is given 
by Bryan (3). The detailed derivations of these equations are adequately 
documented in the cited references and are not given here. | — 
i Solutions t to these buckling equations ¢ can be given as a relationship between 
‘the buckling coefficient k, and the aspect Tatio >,. The subscript b denotes 
that the buckling coefficient and aspect ratio are expressed as functions of 
the of the assembly. the coefficient could 


‘ 


1.—Thin- Walled ‘Members: 8: (a) Stiften Stiffened E Element; (b) intermediately ‘Stiffened 

Element 

ie be expressed as a function of the flat width w, in which w is the width 


_ between the intermediate stiffener and the supported longitudinal edge, so that — 


_ For the local plate buckling mode, the ‘Bryan equation gi gives a ‘minimum ke 


The relationship k, we o, is shown in Fig. 3 for both stiffener 
and local plate buckling modes. These curves are Be similar to those for 


the edge stiffened element shown in Fig.30f Ref.5. 


‘The plate buckling stress coefficients obtained from the buckling equations 
are plotted in Fig. 4 against a nondimensional stiffener parameter A ,/A*. The 


j from Eq. 2 the minimum local plate buckling stress coefficient will be aa 
| 


term A, is s the actual cross- s-asctionsh area of the stiff ener, r, and A? is the stifi fener 

area which would lead to simultaneous stiffener and local plate ‘buckling. GiIajOr 


EXPERIMENTAL INVESTIGATION oti. te 


4 _ The influence of the longitudinal stiffener rigidity on the parsegen pale of 


_ The investigation involved four beam test series with w/t ratios of 47, 70, 
97, and 156, in which w is the flat width between the intermediate longitudinal | 
stiffener and the supported edge of the compression flange. The specimens 


flange, , as shown i in Fig. 5. 


‘The stiffeners of each test series were that moment 


4 
7 


e Buckli jing ng Mode 
FIG. 2. —Critical Buckling Modes 


of inertia of the stiffener for some of the specimens was above and for others 


below the minimum value currently required by the American Iron and Steel 
Institute (AISI) Specification (1). The specific stiffener centroidal moments of 
inertia are given in Table | along with ‘the : specimen dimensions. Table 2 gives _ 
the material properties, experimental ultimate moments, and buckling stress | 

The test setup is shown schematically in Fig-6. 

Strain gages were placed in the central region of the beams to establish the — 

quae stress distribution across the section and ‘to determine when Ic local 
stiffener buckling occurred, or both. +i 


‘The experimental buckling stress coefficients were obtained the strain 
‘reversal method (8) and are compared to those obtained from the ig nl 
_ analysis in Fig. 7. Buckling coefficients for test series I-47 are not shown in 
this figure since the compression flanges of these specimens either did not 
buckle or buckled inelastically. All but one experimental buckling coefficient, 


; shown i in the figure, exceeded the theoretical value for the local eras mode. 


— THIN-WALLED MEMBERS — 
—  Stiffener Buckling 
| 


plausible reason. this is that the compression were restrained 


rotationally by the webs. This was not accounted for in the buckling rer 
-_ _ Nevertheless, agreement between experiment and theory is satisfactory. 
By comparing the longitudinal strain at the stiffener-plate juncture with the 
strain at the web-plate juncture, the relative support provided by different stiff ener . 
_ Tigidities could be assessed. The variation of strain at these locations is shown | 
in Fig. 8 for a representative test series. It is apparent that the strain variations 
at the longitudinally supported edge and at the longitudinal stiffener are virtually 
the same for specimens with sufficiently large stiffeners. The strain variations 
of test I-97-313 (97 is the w/t ratio, and 313 is the nondimensional stiffener 


= of inertia J,/t*) are representative of this type of behavior. In contrast, _ 


gnilaoue 


3 —Critical Buckling Coefficient Curves word 


the strain at the stiffener a specimen with a relatively small 


by test 1-97-115, indicates a stiffener that buckles outwardly in direction 


perpendicular to the plane of the plate. 


of elements that have stiffeners of adequate and inadequate rigidity. In a bi 
a pv section, » formulas will be developed for establishing required stiffener 


of 
af 
= 
] 


‘ for predicting ultimate strengths and for assessing anaes stiffener rigidities 
_ An intermediate stiffener is defined as adequate if the intermediately ‘stiffened 
plate can be treated as an assembly of two stiffened plate elements of width 
w, and a fully effective stiffener. When partially stiffened by an inadequate - 
_ intermediate stiffener, the assembly can be treated as a partially stiffened plate 
4 of width 5 and a reduced or partially effective stiffener. Here b is the width . 
between the ‘supported edges of the intermediately stiffened plate and wis 
the width of each component plate. In the following development, it will be 
assumed that b is twice the width w. This is a valid assumption provided the ‘ 


width of the stiffeneris small relativetob. 


The ultimate strengths of the component plates will be predicted by an effective 
_ width approach. This is consistent with the design procedure developed in Ref. Es 
Sf for edge-stiffened elements. The effective width equation in the current AISI 4 


FIG. 4.—Minimum Critical Buckling Coefficients 
4 specification (@) will be employed. The following i: is a generalization of th 


= 


The effective width of the to total plate of — b can also be > be expressed as 


= 0.95¢ 1 — 0.209 (4) 
Depending upon the degree of support, the intermediately stiffened element 
is categorized as adequately stiffened, partially stiffened, or stiffened. The latter 
is the limiting case, where there is no intermediate stiffener and the plate is 
supported by two webs on its longitudinal edges. These categories are defined 
later, and the assumptions for buckling stress coefficients of each are given. _ 
 Adequately Element.—An intermediately stiffened element is 


iz 
a 
) 


‘@ that the ultimate strength of each component plate element (of width w) 

equals that of an identical plate stiffened by webs on both longitudinal edges. 
‘The effective width of each component plate is predicted by Eq. 3 with 
po a buckling stress coefficient k, equal to four. The intermediate stiffener is 1 
- assumed to remain straight or unbuckled up to the ultimate failure load. Thus, iS 

_ it is assumed to remain fully effective for Purposes of paeeny its contribution ) 

_ Alternatively, the ultimate strength of the total plate of width b can be predicted ie 

Substituting Eqs. ‘San and into Eq. 


i= which i is equivalent to the of ‘each component 


—u 
metal screws ot 


. Partially Stiffened Element.—An intermediately stiffened compression element 
_ is partially stiffened if the stiffener is inadequate; namely, if the flexural rigidity — 
- of the stiffener is insufficient to permit the ultimate strength of each component 


to that of an identical mae by webs on both 


Intuitively, the ultimate ” a stiffened is bounded 
iy two extreme cases: (1) A stiffened plate of width b with an adequate - 
i - intermediate stiffener; and (2) a plate of similar dimensions but without a 
; ue intermediate stiffener. Likewise, ‘the buckling stress coefficient of a ork 


rtially stiffened elements are considered later and will ne Ht 


stiffened element (ks) ps is bounded by the buckling : stress coefficient of an A 
_ adequately stiffened element hous and by the coefficient of an element without | 
an intermediate stiffener (k,),,. The latter two buckling coefficients were shown 
to be 16 and 4 when expressed as functions of the total plate width, " 
| Theoretical values for (k, )ps Tequire a detailed critical buckling analysis. 4 
_ However, for design purposes, (ks) ps may be. determined from the following 
_ expression, which is a close fit (Fig. 9) of "the theoretical buckling stress ‘i 


obtained from a critical bucklin ] 1 
i Titi uckling analysis. A ey 


as b/ ns 


(ks) pe = Fas, thea (8) 


which (ks) os = predicted buckling stress 


Tongs. TABLE 1 —Dimensions of intermediately Stiffened 


in 
inches 
(12) 


“Refer to Figs. 5 and 6. 


elements; J, = moment of inertia of the stiffener. about its axis parallel 
to the plate ‘clement; I* = stiffener centroidal moment of inertia at which critical 
buckling of the assembly is initiated by simultaneous stiffener phe local plate 
buckling; (k,),, = 16; and (k,),, = 4. Eq. 8 was obtained by fitting a quadratic — > 
expression through the theoretical values for (k,),, at J,/1* equal to zero and : 
one. The term /* is the adequate stiffener rigidity based on a criterion of critical fh 
buckling; however, when the flat width ratio of the assembly exceeds some ‘ 
limiting value | (to be determined in a following s ection), adequacy s should be- 


= on an ultimate strength criterion. Since the adequate stiffener rigidity — 
the latter criterion, (J, can be considerably larger than these 


= 
q 
= web, | D 7 
Test b, | in | in in | Lin in | tin | ty in |r, in| dain 
series w/t |inches jinches inches inches |inches| inches |inches inches inches | 
1-47 | 48.0] 6.10 | 5.00 | 9.35] 1.50 | 88. 24. | 0.0565} 0.0575 | 0.210] 0.145 | 60.0) 16.3 
Ta 47.0 | 6.11 | 5.00 | 9.33] 1.55 88. | 24. | 0.0575] 0.0580 | 0.213] 0.223 | 110.0} 19.4 
6.10 | 5.00 9.35 1.55 88. | 24. | 0.0570} 0.0588 | 0.217} 0.259 | 142. | 20.8 J 
47.0} 6.10 | 5.00 | 9.38] 1.54 | 88. | 24. | 0.0578] 0.0572 | 0.209] 0.350 | 217.| 23.5 
46.8 | 6.00 | 5.00 938] 1.50 | 88.| 24. | 0.0577] 0.0590 | 0.215| 0.535 | 499.| 30.2 
«| «71.8 8.89 | 4.08 | 12.0 | 1.50 | 88. 24. | 0.0564] 0.0585 | 0.192] 0.108 40s} 14s 
70.4 4.09 | 120 | 1.50 | 88. | 24. | 0.0571] 0.0585] 0.213] 0.230] 118.) 1998 
70.1 | 8.80 | 4.00 | 12.0 | 1.50 | 88. | 24. |0.0571) 0.0593 | 0.213] 0.230] 118.) 19.8 
69.5] 8.95 | 4.10 | 12.0 | 1.30 | 112. | 37.3 | 0.0590] 0.0593 | 0.217] 0.367 | 232. 24.0 
69.5 | 8.95 | 4.08 | 12.0 | 1.31 | 112. 37.3 0.0592 | 0.224] 0.458 | 367. | 27.6 
9.8} 8.95 | 4.10 | 12.0 | 1.34 | 112. | 37.3 | 0.0587] 0.0592 | 0.224] 0.532 | 488. | 30.1 
69.3} 8.93 | 4.00 | 12.0 | 1.35 | 112. | 37.3 | 0.0590] 0.0598 | 0.230] 1.150 |2650.| 51.2 
| (96.4 5.50 | 15.0 | 1.56 88. | 24. | 0.0570] 0.0690] 0.219] 0.226 | 115. | 19.6 
11.8 | 5.45 | 15.0 1.58 88. | 24. | 0.0564] 0.0690] 0.214] 0.263 | 150.| 21.2 
79.3] 11.7 | 5.50] 15.0 | 1.30 | 24. | 0.0694) 0.0698 | 0.236] 0.359 | 153. | 21.0 
a 79.5 5.50 | 15.0 | 1.50 88. 24. | 0.0692] 0.0697} 0.224] 0.545 | 25.9 | 
80.7 | 11.7 | 5.50 | 15.0 | 1.55 88. | 24. | 0.0681) 0.0695 | 0.228} 0.681 | 527.] 
185. | 186 | 5.50 219 1.55 88. | 14. | 0.0580} 0.0595 | 0.217] 0.218 | 107. 19.3 
156. | 18.6 | 5.49 | 21.9 1.50 | 14. | 0.0580] 0.0573] 0.218] 0.268 | 147. 21.0 
159. | 18.7 | 5.50 | 21.8 | 1.54 88. 14. | 0.0572] 0.0575 | 0.218] 0.350 | 235. a 
157. | 18.6 | 5.50 | 21.8 | 1.57 | 88. 14. | 0.0573) 0.0575 | 0.217 0.514 477. 
&§ 


1 


= 


criteria are not ide identical. Consequently, it would not be appropriate to replace — 

For this reason, Eq. 8 is restated in terms of (J, Denetiisdy and an exponent > 
a n, in which n will be chosen such that the following oe is a close a b 
approximation of the theoretical buckling coefficients: 


_ The effective width of a partially stiffened ed element of of width bis then d determined 


Material Properties 


Results 
inkips per | M,,,in 
square inch _ kip- 

(5) 


sur 
G 


. 


“Inelastic buckling stress coefficient. = 
No observed stiffener or local plate buckling. 


by Eq. 4 with (k,),, given by Eq. 9. The contribution of the partial intermedia ie 
7 stiffener to the ultimate strength of the assembly is covered in the following — 


Bs _ The procedure for predicting the contribution of a partial or buckled stiffener 
= Ge ultimate ultimate strength o of the s assembly, i ist the same as outlined for edge- stiffened 


| | 
: 
/ ‘Test in kips pe 
= | | | we] 
40 | | 159 | 16.7 
107. 2.40 | 509 | 104. i: 10.6 
425 | 30.2 | 106. 
=: oo | 504 } was | m2 


elements in Ref. 5. The cross- -sectional area ea of the stiffener is reduced — 
by a simple linear expression to reflect its reduced sateetetimaisiath in in resisting — 


in in which 4 the unreduced cross-sectional area of the stiffener; et” 
moment interia; and (J,) = the minimum 
necessary to support the. plate adequately. 
_ Experimental justification for Eq. 10 is given in n Fig. 10, weal the experiniaiel 
data is shown to scatter about the linear relationship given by the foregoing 
equation. The experimental partial stiffener areas for this figure were obtained — 
by subtracting the contribution of the predicted effective areas of the webs _ 
and from. the total effective section modulus the beam specimen at 


FIG. 


In limit, the ‘partially stiffened element adequately stiffened, 


In this a. the effective widths for adequately and partially stiffened Gonaete 


_ §$tiffened Element.—A stiffened element can be considered as the limiting 


case of an intermediately stiffened element. Therefore, the equation for predicting _ 
_ effective widths of elements partially stiffened by an intermediate stiffener must 
also be consistent with the effective width approach currently employed to 


this case, equals zerc zero and 9 reduces t to 


t which ( ok , = 4. With sian ie for k,, Eq. 4 becomes the effective width ; 
_ equation used in the AISI specification for stiffened elements of widthb. 


a 
| 
] 
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_ The definition of an adequate stiffener could be based on either a critical | hes 
_ buckling criterion (CBC) or an ultimate strength criterion (USC). With the CBC, al 
the adequate : stiffener is defined as that "stiffener rigidity at which instability _ 
3 initiated by simultaneous stiffener buckling and local plate buckling. With 
the USC, the adequate stiffener rigidity is defined as the minimum rigidity 
g at which the ultimate strength of each component plate on either side of the 7 
stiffener equals that of an identical plate stiffened by webs on both edges. os 
a u CBC. —To obtain a stiffener requirement based on a CBC, the stiffener Sing vi 
equation is solved for those stiffener dimensions for which the minimum buckling : a 
stress coefficient k, equals 16. This is the stiffener dimension for which eae 
buckling and local plate buckling | occur simultaneously for all aspect ratios , 
The requirement is shown in Fig. 11 as (J,/t* hee the required stiffener rigidity 4 
nondimensionalized wah Respect to the plate thickness ¢. Three 


= 


Stiffener Buckling Mode Series 


4 


are shown in “this figure. For 1 the larger requirement, it is | assumed that - 


neutral axis of bending of the stiffener coincides with its centroidal axis, and 
for the smaller it is assumed that the neutral axis of the stiffener coincides 
with the middle surface of the plate. | Neither of these, _ however, is ‘strictly 
correct. The stiffener attempts to bend about its centroidal axis, and the plate 
attempts — to bend about its middle surface, as though each would buckle 
independently of the other. Because of compatibility of strains, , shear stresses 
Se at the stiffener-plate juncture and the actual stiffener neutral axis hog 


- peflects this change in neutral axis. The third stiffener requirement, shown in 
: Fig. 11, is based on a stiffener moment of inertia due to the adjusted neutral 

axis. The stiffener requirement based on an adjusted neutral axis is not — 
Significantly different from that which assumes that the neutral axis is located 
at the plate’s middle surface. ‘Some design specifications (1) implicitly make 
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= assumption and for ae case it appears to be valid. Nevertheless, the CBC 
stiffener requirement for intermediate stiffeners (and for edge stiffeners in Ref. — 

will be based on the adjusted neutral axis; 
_ To assess the applicability of the CBC requirement, the observed ultimate 
"strengths of several series | of longitudinally stiffened elements are “compared 
in Fig. 12. Each series consists of specimens with identical “plate dimensions — 
e different stiffener rigidities. An increase in stiffener rigidity beyond that 
¥ required by the CBC increases the ultimate strength of the assembly to varying — 
For test series ‘1-47, 70, 1-97, and I-156, they are about 9%, 15%, 
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FIG. 8.—Post- “Buckling Behavior of Stiffened Assembly 
mate? 


ae to some increase in ultimate strength; however, this increase amounts 
_ toa few percent at most. The significant i increase in ultimate a is attributed © 


es Also, it is s apparent ‘from Fig. 12 that for | test s series s 1-47 ‘only small increases f 

in ultimate strength are observed when the stiffener rigidity i is increased beyond — 
the cee requirement. For such small w/t t a, ‘this implies that 


G 
22. 
j 
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as stiffener based of on the CBC is ‘not ‘significantly different peng 
a. based on the USC; therefore, only for w/t ratios in the advanced post-buck- <a 
- ling range is the CBC stiffener requirement excessively unconservative. The 
_ | a In the range of w/t ratios for which the component plate of width w remains — 
i effective (hereafter called the fully effective range), the stiffener rigidity 


based on the CBC requirement provides the minimum required support so that 
stiffener buckling occurs simultaneously with material yielding. Consequently, 
= @€—s the stiffener remains unbuckled up to the ultimate failure load. The CBC and 
stiffener requirements are identical in this range. 
On the other hand, for w/t ratios larger than the limiting w/t ratio (the 
post-buckling g range) the stiffener tigidity based on the CBC requirement provides — 
support such that simultaneous stiffener and local plate buckling occur before 
material yielding. For this range, the stiffener will deflect normal to the plane E ae 
oft the plate prior to the stiffener-plate assembly reaching its ultimate load. a 
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9. —~Veriation of Stitfener Ri Rigidity 


Consequently, the _— requirements based on the CBC and USC will dif fer 
Inthe following section, stiffener requirements based primarily on experimental 
7 results will be developed for the usc. It will be ‘shown that for w/t Tatios 


is insufficient for purposes of predicting ultimate strength. 
al USC.—In Ref. 5, the adequate edge stiffener requirement \ was expressed as 
- a function of the w/t ratio normalized with respect to (w/t). in which (w/t), 


from Eq. 3(b): sae 


| 

IS Lhe Tatio DeIOW WHICH Lhe Ol Width W 1S CLICCLIVE aS an 


 forsteel with k,, equal to 4. By no normalizing the stiffener 
F with r respect to (w/t), , the requirement can be expressed as one equation rather _ 
than a family of equations, one for each yield stress, iki 


— 


= For edge-stiffened elements, the transition between the fully « effective a) = 


and the p post- -buckling was yas taken as 


_ For the range w/t: < (w/t), the critical buckling criterion was used to determine 


IG. 10-—Eetive Stiffener Cross-Section 


Therefore, ‘ee transition between the range of the CBC and the USC will be 
a. extended to (w/t)/(w/t), = 1.5 instead of 1.0, and the required minimum moment — 
of inertia can be adjusted to agree with the adequate rigidity based on = 
e series I-47. For lack of experimental results in the range of (w/t)/(w/t), ratios — 
i. below I. 4, the ‘stiffener requirement i is varied linearly in this 1 range. ge 7 


—— the value of the adequate stiffener rigidity. For the intermediate stiffener 
requirement, a slight modification of this is warranted. For the test series 1-47, 
h 
| 


ik 


‘is which (w/t), i is defined b = Eq. 13, and d (w/t) i is the flat width re ratio below 


11.—Stiffener Requirement Based on 


221, 


Eq. 18 can mbes ‘expressed in terms of 9 if ubstitutes b = 2 w. yy se 


For the range w/t > 1.5 (w/t), , aprocedure similar to that for the edge-stiffened 4 ; 
elements (5) is employed. Eight hypothetical stiffener requirements, A—H (in pe 
_ Fig. 13), will be compared to the test results. These requirements are arbitrarily 


drawn straight lines of i increasing slope havin — a common at at(w w/ /(w w/ 


(2) 
ke adequate 
i 
ae J Neutral Axis at Centroid 
J 
q 
¢ 


sts THIN-WALLED MEMBERS 


- equals 1.5. The requirement for which the protien and enpedenuanal ultimate 4 


strengths are in best agreement will be suggested as a possible design requirement. 

oF Previcteo AND ExPeriMENTAL Uttimate STRENGTHS 


To facilitate a ‘comparison mn of predicted and experimental | ultimate strengths, 
two nondimensional ratios are used oF 


aw which (Padi Predicted ultimate by the 


a. — FIG. 12- —Variation of Ultimate Strength Versus Stiffener Rigidity 


outlined for adequately or ‘partially stiffened elements (whichever i is uote 
(Px), = experimental ultimate strength; and (P,,,),, and (P,,),, = predicted 
ultimate strengths for an adequately stiffened element and an element with 
no intermediate stiffener. All predicted ultimate baroy are determined by 
the effective width approach previously outlined. 
These ratios normalize the predicted and ultimate strengths such 
that they provide a measure of the degree to which the intermediate stiffener 
4 supports the flange. That is, if the ratios equal one, the ultimate strength of 
the intermediately stiffened flange equals that of an adequately stiffened flange; 
an and if they equal zero, » the ultimate strength equals that of a stiffened element 


641 
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without an In the following, experimental and p 
wimete strengths are compared using the previously defined ratios. _ sedate 
Range w/t, < w/t = 1.5 (w/t), .—For this range only one test series was 
> conducted. For this series the ratios of experimental to predicted ultimate strengths — 
R,/R are given in Table 3 and are graphically compared in Fig. 14. In this 
= the adequate stiffener requirement is given by Eq. 16, and (k b)ps IS 
determined from Eq.9withmequals2, ows 
- Correlation between experimental and predicted ultimate strengths is satisfac- 
tory, which is to be expected since the required stiffener rigidity was adjusted 
based on this test series. It is noted, however, that the predicted ultimate strength © 


: two adequately stiffened tests are overly conservative. This is s attributed 
* the partial plastification of the webs of these beam specimens. It has been 
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FIG. 13.—Hypothetical Stiffener Requirements: usc 
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smaller than, (w/t), have ar an inelastic r reserve * capacity duet to: a pastiol plastification 

_ of the webs. Since this was not accounted for in the analysis, it is in all likelihood 

the reason for the predicted conservative ultimate strengths for several of these 

; Range w/t > 1.5 (w/t), _.—For this range, (k,),, is given by Eq. 9 as a 


function of (J, wont and n. In the following comparisons of experimental il 
predicted ultimate strengths, eight hypothetical stiffener requirements are as- aly 
sumed for (J,).sequare and several values for m (1-5) are considered. The | 
comparisons of R, and R have been extended to determine both a andl $n 


stiffener requirement and a —- approximation of of (ks) os . for that for that requirement. ba 


virement and a g 


4 
= 


=o ratios R,/R are compared. To identify my fe value of n in 1 these comparisons, 

gE Ri is subscripted by the appropriate integer n. Also, to determine the combination 
«(Of coefficient which best geodicts the experi- 


mental Ultimate Strengths is in Range (w/t), = 5 
ars 
4 
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Test Series 1-47 


FIG. of Experimental and Predicted Ultimate Strengths: < 


mental ultimate strengths, it is arbitrarily assumed that a predicted ultim 
_ strength agrees satisfactorily with the experimental i ee 


20 
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TABLE ee _Number of Ultimate Strengths Accurately Predict 


6 Re 
— Stiffener 


requirement 


bale 


a : 
“The largest number of accurately ultimate strengths. 
Note: 


"TABLE 5. —Comparison of Ultimate Strengths i in pies w/t = 1.5 (w/t),: Intermediate 


= Stiffened: 


adequate 


Test | | 


series 


1.53 
4.57 


T4 
ed: Range w/t = 1.5. 
ie 
Me, 3 
w/w, q R,/R, Ro/ Ry R,/R; R,/R 
1-70 | 70} 40.5] 0.151 |2.10 |0.500] 2.22 [1.18 |0.952 |0.855 0.806 
4 18. 0.480 | | 0.682 | 1.22 |0.990 |0.926 0.893 [0.877 | 
232. 0.937 | | 0.887 | 0.935 |0.917 |0.917 0.909 0.909 | | 
97. 115. | 0.256 |2.87 | 0.593} 1.77 |1.00 j0.893 
80. 153. | 0.461 | 2.37 [0.704 1.29 |1.04 0.971 |0.935 |0.917 
| 0.956 ons 0.877 |0.870 0.862 
| 156. 0.500 2.73 [1.35 |1.06 0.935 |o.870 
| | 0.486 2.10 |1.16 |0.952 | 
| 0.256 | 10.563} 1.70 1.12 j0971 Jo901 | 
| 


The: Situ of tests that : satisfies this inequality i is given in Table 4 for 5 od 
combination of stiffener requirement andm, 
Po ‘From this comparison, the experimental ultimate strengths are best het 


i - requirement F, since the largest number of tests satisfies | the foregoing al 
lit 


ity for this requirement. The requirement is 


Ro Rrat ration are for ener F. The arithmetic 
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20 40 BO” 
adequate 


02 
Requirement 


2040 80160 02 04 06 08 1: 

15.- —Comparivon of of Experimental and Predicted = 
Hee. are 0.961 and 0.05. For partially stiffened elements, n = 3 is the best choice << 
1 2 for n, since the predicted ultimate strengths are pre Sea conservative for — 
= 3 compared to n = 4. The arithmetic mean and standard deviation for ; 
are 1.02 and 0.13, which are quite satisfactory. 
The and predicted R ratios are compared for all 


¥ 
a 
- 
| 
ae 


1981 


the stiffener. requirements in Fig. 15, which confirms that the ¢ ain requirement 


-F is the most satisfactory of the eight. In each plot of Fig. 1S, the upper + 
curve i is for w/t = 70 and the lower curve is for w/t = = 156. ¥ = 


‘The approach developed here is limited to intermediately stiffened a 
_whose stiffener width is small relative to the width of the assembly (i.e., 

approx equals 5/2 in Fig. 5). When w is smaller than 5/2 
equations for adequate stiffener rigidity and effec tive may be inappropriate, 
_ Two reasons are given for this: 


® 


04 06 os 10 20 02 a4 08 
FIG. 15.—Continued 
= * The critical —_— equation given in Ref. 2 for partially : stiffened elements 
(and on which the aforementioned development is in part based) i is not ae 
_ to this case. In that derivation, it is assumed that w equals , ar se 
the behavior of the assembly. 
These influences are currently under study at Cornell University with the objective ~~ 
extending the approach developed here. ite wat shania. 


SumMaRY AND Conctusions 


«TSA 


"adequate 


tina 
_ An approach has been presented for predicting effective widths of interme- 
diately stiffened compression elements that are either t adequately or partially 


1 

= 
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stiffened. / A stiffener on minimum required stiffe 


rigidity to support these elements adequately is also presented. Bi es 
Two criteria were considered for the stiffener adequacy. For flanges 
_ which are in the range (w/t), < w/t = 1.5 (w/t),, a stiffener nelenes. 
_ based on a CBC appears satisfactory. For ated having larger w w/t t ratios, 
a requirement based on an an USC i is } mecessary. 
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are used in this paper: 


‘tee ny cross-sectional area of intermediate stiffener; 


The a mis = characteristic stiffener cross-sectional area at which stiffener 
buckling and plate buckling occur simultaneously; 
effective cross-sectional area of partial intermediate stiffener; a 
= total width of intermediately stiffened assembly; aoe . 
= effective width of plate components of intermediately stiffened 


_ limiting width of stiffened plate with adequate size inter 
stiffener below which it is fully effective; 
limiting width ofa simple stiffened plate (i.e., 


with no intermediate a aiaatedh below which it is fully effective; 


rican 
g | 
— 
Mass 
&§ ol lement 
— 


Fe Young’ s s modulus = = 30,000 ksi 4 210,00 000 ) NM/m 
moment of inertia of edge stiffener taken about | its centroidal 
* = characteristic stiffener centroidal moment of inertia at which - 
¥ ocal plate buckling of flange and stiffener buckling occur 
a | plate buckling of flange and stiffener buckling 
pore? 
required minimum stiffener moment of inertia to 
support flange; 
buckling coefficient for plate of width 5; 
buckling coefficient of intermediately stiffened assembly that 
= : buckling « coefficient of intermediately stiffened assembly that — 
is partially stiffened; 
buckling coefficient of simple stiffened ‘element of width ‘b; sal 
= buckling coefficient for plateof widthw; 
L- supported length of beam specimens; = 
from concentrated ioad to support in beam tests; au 
pp 
= experimental ultimate bending moment of beam specimens; r 
Oda "predicted ultimate strength of intermediately stiffened assembly 
with adequate size stiffener; 
= predicted ultimate strength of "adequately, partially, or or simple 
stiffened element (whichever i is appropriate); 
was = predicted ultimate strength of simple stiffened element 
CO stiffened element with no intermediate stiffener); 


(P experimental ultimate strength; 


normalized predicted ultimate strength; 
= normalized experimental ultimate strength; ” 
’ s width of component plate of intermediately stiffened assembly; 


effective width of plateof widthw; 


- limiting width of flange below which i it is fully effective as 
stiffened element of width w; fom 

=k E/12(1 — (w/t)? = 


= material yield stress, in kips per square inch ksi : =6. 9 
flange aspect ratio. 
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-PRoposeD STEEL COLUMN STRENGTH rH CRITER RIA 


For 200 the si strength of metal columns has 1 attention. 
The Euler formula derived for elastic columns has been shown by tests to 
_ give an upper bound for rather slender columns. The inelastic range of column 7 
: slenderness has proven much more elusive to quantify. Strength of such columns ~ 
- depends strongly on at least two parameters, residual stresses and crookedness, 7 7 
not quantifiable by the designer. Accordingly, it is not surprising that the scatter 
_ in the inelastic range of column strength is larger than that for most other “yy 
It was not until the 1940s that Shanley explained the inelastic column failure 
_ ‘mechanism as we presently understand it (19). . This mathematical solution | for 
column strength i is iterative and best handled by computer (1,2). Shanley’ s model | 
_has been used to accurately predict the strength of many column tests if sufficient 
_ detailed information is provided. Because of the lack of a closed solution at 
this time and the large scatter in strength, column design equations must be 
empirical in the inelastic range of slenderness. The basis of the equations , | 
either physical tests or theoretically calculated column strengths. 
us _ This paper presents the development of empirical column strength criteria 
a based on physical tests of steel columns. These tests are the only sound evidence 
that includes all of the variables in real steel columns. However, most of the 
_ tests have pinned ends; thus, the tests do not simulate an actual column in 
‘The test data can be fitted with the Euler hyperbola in the elastic ra range. = 
a The inelastic range of column strength can be represented by a sloping straight _ 
line from the Euler curve to a horizontal plateau in the lowest slenderness. 
zz column data are further examined to provide a compression resistance 6 
factor, ®., , which can be used in the load and resistance factor design a 


pr. Consulting Engr., Technical Services Div., Engrg. Dept., B Bethlehem Steel Corp., 
-Note.—Discussion open until September 1, 1981. To extend the ¢ closing date one month, < a 
a written request must be filed with the Manager of Technical and Professional Publications, — 
= ASCE. Manuscript was submitted for review for possible publication on August 13, , 1980. 
This paper is part of the Journal of the Structural Division, ye of the American © oy 
Society of Civil 1981. 0044- 


i 
| 
| 
| 
| il 
| 
| 
: 


APRIL 


i. of steel buildings. The resistance factor is a measure of the reliability ol 
assigned to the nominal strength The smaller scatter o 


ults, to .0 is the. resistar 


The n most common of a cross- configuration 
_ which generally buckles in the pure flexural mode. These sections include ‘tolled 
7 ‘I-shapes, rectangular boxes, and solid and hollow rounds. Other shapes, such 


‘ss angles, are subject to flexural-torsional buckling and are not | considered. 
Also, sections which have critical local buckling of elements of the cross section 


are not considered inthis study. 


The lowest slenderness columns fail by yielding in compression. The strength © 
a these stub sections is a function of cross-sectional area, compressive yield 
- stress, and strain hardening. Columns in the intermediate range of slenderness 

_ have strengths dependent on a number of additional variables, including slen- : 
derness, initial tensile yield stress, elastic modulus, moment of 


a Flexural buckling of columns can be divided into three regions of slenderness. 


] 


inertia, residual stresses, , and end conditions. Columns with large slenderness ; 


‘The test data used in 
-_ageregations of tests performed over the past 30 3 yr by various investigators. 

of these data were cataloged and reported in Ref.12. | 


Column test t data may be examined be plotting average stress 


columns having different yield ‘it is common nondimensionalize the 
4 strength and slenderness by the yield stress. The compressive stress at failure — 
is divided by the yield stress. Slenderness is transferred to d » in which = 


V F, E). The transformation of slenderness i is derived from the elastic 
_ buckling equation. In this study, the test data were nondimensionalized in this 
manner using laboratory-measured tensile yield stresses. 
Table 1 gives a breakdown, by column type, of the 393 tests included in the | 
_ data base. The rolled H-shapes in parentheses were part of a single major 

test effort. The scatter in these data was less than in the remaining data, aden ; 
that test methods may have introduced some of the observed scatter. 
The data were divided into two categories—columns with //r less than 110, 
and those with slenderness of 110 or greater. The more slender columns were 
nel to behave as elastic members and fit a hyperbola. A plot of the data 
t Fig. | indicates this to be true. The less slender columns exhibited greater 


scatter, and the best curve type was not apparent. 


In order to select an appropriate curve for the inelastic range, test data of 

shapes with the maximum residual stress reported were selected from the 393 
4 Forty-four such column tests were identified. Appendix I shows a summary 
of a statistical analysis of 15 various commento of four parameters: nondimen- 


sional slenderness, measure yield stress, ‘maximum measured compecssive — 
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PROPOSED | STRENGTH CRITERIA 
=f _re Stress, and shape of the cross section. Based on these results, a a 


straight line relationship between strength and slenderness was 


selected because it was statistically as good as the 


point of view it was the simplest. ya 

TABLE 1.—Column Data Base ube 


Flame cut 
As-roliled 
Cold straightened 


d _ 


=i 


be “22 tests were plotted; information taken from Refs. 4, 6, 13, 14, 17, 21, and 22. 


tests were plotted; information taken from Refs. 17 and 26. 14 
“14 tests were plotted; information taken from Refs. 4, 14, and 26. Ae, a ni 
*22 tests were plotted; information taken from Refs. 1,9 15, and 25. 
tests were plotted; information taken from Ref. 16. 
£393 tests were plotted. 


Note: Ref. 12 was used as as the asic : eference to build this data base. 7 


~~ : column tests with nominal yield stress between 33 ksi and 36 ksi (227.7 
* a. mPa and 248.4 mPa), and column tests with nominal yield stress of 100 ksi 


: 
- p 
Totals? 
tm 
Notapplicable | | | | 4] 2 | 2% 
Se 
| 
| 
q 


= 


AP 


(690 mPa). No columns were tested with slenderness greater than 110 for 100 - 


‘This division by yield stress ¥ was made to investigate the ‘possible effect. of 


- yield stress not accounted for correctly in the nondimensionalization process. 


Maximum Stress/Lab. Yield Stress 


tensile yield stress was weed to each test 
Table 2 gives pertinent statistics from the regression int The followin 


Ff = 33-36 ksi (227.7- 2468. 4 mPa) 


06 (068 10 12 14 «16 24 26 


Po , with a coefficient of determination, R’, of 0.74, explains 14% of the — 

observed variation, which is better than any of the equations reported in Appendix _ 
‘I. One reason for the improvement is that Eq. | is based on a number of 
~~. with many replications, whereas the equations in n Appendix 1 were based 


errors were associated with test specimens of low yield stress. Some. of these be 
specimens had reported laboratory yield stress as low as 28 ksi. This is not 


| 
0.2 
0.204 
=: | 
| he or the three fitted 
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PROPOSED STRENGTH CRITERIA 


Eq. 2 is a hyperbola v very similar to ‘the Euler equation. on. The coefficient of , 
determination is 0.96. This increase over the R* for the inelastic region is to 
t expected because elastic buckling, which is dominant in this region of 
ethereal fewer parameters affecting strength. Also, calculation of \ 
is theoretically correct in the elastic region, whereas, the use of lambda . 


7 the inelastic region may introduce scatter. Eqs. 1 and 2 are plotted with the 


data used in their development in Fig. 1. 

% Eq. 3 is for columns of material with nominal yield stress ; equal to 100 ksi 

: (690 mPa). There were 36 data points in this category. The coefficient of 
determination is 0.77, which is similar to that for Eq. However, the largest 
— error sl been reduced to ) 23%, from 45% for Eq. 1. This reduction is, in 


_| Square inch coefficient 


2 


part, due to the smaller scatter erin yield ield stress. | Eq. 3 is p ‘the 
used in its development in Fig.2, “4 


m, The similarity between Eqs. 1 and 3 leads to the. decision to combine columns 
of nominal yield stress material from 33 ksi-100 ksi (227.7 mPa-690 mPa) into © 
single combination results in Eq. 4. 


q 


_ columns of both levels of yield stress 33 ksi-36 ksi and 100 ksi (227.7 mPa-248.4 
mPa and 690 mPa) i in the inelastic range of slenderness. Eq. 4i is tangent to 
the Euler hyperbola at \_ =1. 53. It has been modified to account for a 7% — 
difference between actual average yield stress and minimum specified yield ; 
- stress. This modification is shown in Appendix III. For slenderness greater 
a _ ‘than 1.53, the Euler hyperbola represents column strength. The fitted hyperbola 
: in the elastic range is modified to coincide with the Euler curve. The mathematics 
of this modification is shown in Appendix II. For slenderness greater — 


1.53, the Euler hyperbola represents the nominal column strength. x 


A further modification is made in the low slenderness neat because Eq. 
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4 predicts a a of 1.3 times the yield strength Although 
a tests indicate the strength of some columns may be this high, it is doubtful ‘“ 
that such a high value represents typical behavior. Local buckling limits in ol 
the code are based on yield point. _Galambos (10) reports that the mean value 


is lower than the scatter of ‘column strengths in “the low slenderness range. 
_ This difference is likely due to strain hardening and possibly higher compression = 
‘yield stress. A plateau has been placed at 1.10F, for steels with nominal yield 
stress of 65 ksi or less to recognize some benefit for overrun of yield stress ys 
and strain hardening. Steels with nominal yield stress over 65 ksi do not exhibit o 
the -— overrun in yield stress. Therefore, the plateau for these steels has 


o 


he 
Fy = 100 ksi (690 
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wing are nominal column equation 


= 


< 65 ksi; 


035<A< 1.53 


= a shape fact factor which accounts for certain cross section proportions = 


| 
= 1.0QFA,, 052 <0.53 for > 65 ksi 
=|13-057(—)\/ 2 lor ......6 | 
m =| 1.3-0.57(—] | OF; 4,, 


and is equal to 1.0 for on’ shapes rested. Refs. 3 20 one 27 « cover Q. Because 
these equations are intended design use, they are in ‘terms of the 
minimum specified yield stress, F;. The term R,. is the nominal resistance 
ay _ The LRFD model ‘employed by Galambos and Ravindra (10) utilizes the 
coefficient of variation (COV). The COV is determined by dividing the standard 
z deviation by the mean value, o/. In many instances, the variance is constant : tJ 
_ (over the range and the mean differs. The result is that the COV is not constant. 
In the column case, the mean value decreases with increasing slenderness. Thus, 
the COV becomes larger with increasing slenderness if standard deviation is fj 
constant over the range. This problem is evident in Table 3 which presents 
a summary of calculations used to determine a weighted COV for the column 
_ data based on the fitted curve. Calculated mean values of strength at midrange ~ F 


TABLE 3 —Coeticient of Variation from Data i in Inelastic 


: jen 
0.35-0.55 
0.56-0.75 
0. 76-0.95 


“Weighted CW 556/258 = 0.0913. 


\ are used as the mean for each \ segment. The sander’ error of the estimate — 
has been used as an estimate of the standard deviation as done by Galambos. 
Comparison of the weighted COV indicates it is a fair representation of the 
The compressive resistance factor, ®., is calculated in Appendix II for the 
entire range of slenderness. The COV in the elastic range is smaller than that 4 
ag sce in Table 3 for the inelastic pa, As a result, the resistance factor, _ 
, is larger. For the sake of simplicity, the hyperbola has been shifted upward 7 
"coincide ' with the Euler | curve, as ‘previously mentioned, and the resistance 
- factor has been reduced a like amount. The result is a constant resistance factor = 


‘The regression analysis has comt combined all column types. Each type of column 

included in the : set will be compared to the nominal strength curve and this — 


“curve reduced by the resistance factor, ® 


| 
t 
| 1.06-1.25 | 1.10 | 4 0.673 0.087 | 7 6.699 
1.35 | (0.531 0.060) 22 (2.618 
7 | 


_ The results of 324 rolled shape column tests shown on Fig. 3 are compared © 
to a plot of the proposed nominal column strength, and ®. times the nominal 4 
column strength. The data plotted in Fig. 3 include seven tests of W12 x on 
and 16 tests of their European counterpart, HEM 340. The slenderness of oo 
23 columns is either 50 or 95. Their laboratory measured tensile yield stress g 
was between 27 ksi and 35 ksi. The results of all tests of rolled shapes, except — 
five Wi2 x 61 tests at //r = 50, fall above the ®.(R,.) curve. The points 
below the Curve were tests performed on American shapes thet were 


= not have been sold for structural members od 
11, 16, 17, 18, 22, 23,2426 


y-Axis)= 33-3 36 ksi (227.7-248.4 mPa) 
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_ Fig. 4 shows test values of the W12 x 161 sections with //r = 50 mete 
ie maximum initial crookedness. These tests were performed using two 


different loading methods. The earlier tests were performed by centering the 


midheight of the column for low loads. “This method of loading i is called rad 
“‘old Lehigh method.” The more recent tests were performed with the loads 
applied at the geometric center at the ends of the column. The figure dramatically _ 
demonstrates the effect of both initial crookedness and test method on column — 
4q 
strength. Those specimens tested by geometrically centering the load at the | 
; ends of the column are further divided between dynamic and static tests. Test 7 
a a values determined at both testing rates were nondimensionalized by the same _ 


laboratory determined static tensile stress. pats 
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a Figs. and 6 ‘residual ‘stresses in ‘the. American W12 x 161 
section which has not been straightened and the similar European section, HEM | 
340, which may have been straightened and thus is representative of normal 
practice. The HEM 340 sections all had test strengths greater than ®.(R,,.) 
These HEM 340 sections were tested with the load applied to to the ‘geometric 
center of the cross section at each end of the column. 7 _—— 
- Literature provides little specific information concerning the effect of on ¥ 
straightening on the strength of columns, although crookedness has long been — 
identified as a contributor to premature column failure. Cold straightening has 
two beneficial effects. First, the straightened column is subjected to smaller 7 
bending moments due to crookedness. Second, cold straightening tends 
to reduce any longitudinal residual stresses due to mney cooling during 


or or due to cold forming. 


x Static test 


= Average static yield — 


e= Initial mid height eccentricity 
= Maximum test load rw 


Actual gross area of column. 


straightness 
limit .001 
Ref.20 

Dor 


ASTM AG out-of- Wig 


Z/R = 50 


4. 4.—Column ‘Strength V initial Crookedness 
is performed on the majority of steel product less than 100 lbs/ft (149 kg/m) 
; ¥ _ and 18 in. (457 mm) deep, although the capability to straighten different shapes a 

varies between mills. Rotary straighteners bend the section back and forth over 
a sharp radius, causing plastic deformations which obliterate previous residual 
stresses. This is a continuous process affecting the entire length of product. 4 
The residual stresses resulting from this process are reported to be rather small. be 
_ gag press i is used t to ) straighten large sections that cannot be poral 4 


Small sections are sometimes straightened by stretching, which also el eliminates 


| 
| ] 
| 
Ee. 
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‘residual process is is often used on aluminum sections. 
‘The tests of rolled shapes, except the excessively crooked atypical W12 x 
F 161 sections, confirm that the proposed design criterion is safe. The unstraightened 
shapes confirm the rationale of limiting crookedness of steel compression 
members, even in the rather low Tange of s slenderness. oben ie SMaTT = 
Box-Sxares ano H-Suares 


Tests of 22 welded H-shapes and 14 welded box shapes are plotted with — 
the R,. and ®.(R,..) curves on Figs. 7 and 8, respectively. These test data, 
“except | for three points, lie above the ®.(R,..) curve with the central tendency 
somewhat below the R,. curve. This plot confirms that the design criteria are’ 7 
compatible 1 with the physical evidence of welded column shapes, except for 


Reference 23. 


a 
4 


© Inner 


and A7 steel was the common » strecteral grade. The welding w was done a 
using the submerged arc Wire size was 5/64 in.; 


_ input would most likely have caused distortion of the ‘section, resulting in 
subsequent straightening which would have modified residual stress patterns. 
_ The actual specimens were reported to be very straight, with maximum initial 
crookedness less than 0. 0. 001. ‘Modern welded H-shapes would higher 


| 

“4 

ore tresses in W12x 161 

a the three low values for tests of H-shapes which deserve further investigation. - | 

P 
three low points were tests of the shape shown in Fig. 9. These tests 

J 
d _ The relatively low heat input probably deposited considerably less weld, and | 


Although there were three low values for some of the earliest welded columns, — - 
‘the proposed criteria safely represent the strength of modern welded columns 


AND Coverptatep H-Suares 


The test data for nine riveted columns shown in Fig. 10 fall around the nominal — 


strength. The | ‘scatter, although substantial, seems to be less than that of ‘the 
welded columns. ‘This may be due to the smaller variation of residual stresses, 
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FIG. 6. —Residual Stresses in HEM 340 (Italy) 


coverplates reported are also shown in Fig. 10. ‘These types of columns may _ 
also be safely represented by the proposed cr criteria. a 
Tests of 22 solid round bars are presented in Fig. 11. All tests are of 100- ksi 
- (690-mPa) yield stress material. Unlike the other column tests with pinned ends, 
4 these tests were made on flat end specimens. Strain gages were used to determine : 
_ the inflection points so that effective slenderness could be determined. The 
| § effective slenderness was used to plot points on Fig. 11. These points are — 
well-distributed about the fitted line, indicating that the strength of this c column 
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Comparisons oF Resutts Tot Orner 

7 5 Fig. 12 compares pares the ® ie .) curve developed i in this work with the ® oe 


curve proposed by Galambos (8) and the Canadian limit state specification (3). 

~ 
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Additional abscissas of l/r for F, equal 36 ksi and 50 ksi (248. 4 mPa and 345 
mPa) are also shown. The proposal by Galambos results in nearly | a straight — 
line which lies somewhat below this proposal. The Canadian limit state eueden 


crosses the proposed curve several times, but closely parallels it over the entire a 


TI 

P 
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PROPOSED STRENGTH CRITERIA 

- slenderness range. e. The proposal by Galambos, which was based on a theoretical 
model of column behavior rather than test results, ‘indicates that a straight line 
is ar reasonable representation of column strength. Among the reasons why 


ye 
Laboratory Measured Yield Stress = 335 
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pongetegy FIG. 10. _—Coverplated H and Riveted Test Data 
the Galambos curve is lower is the fact that this model was not adjusted for 
the actual yield stress of columns which is reported to be about 7% higher 


than minimum specified | for hot-rolled shapes. A second reason ‘that it is 
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FIG, 9. —Residual Stress Pattern in 
10 + Fy = 33-96 ksi (227.7-248.4 mPa) 
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based on the conservative assumption of all een: aving the ‘maximum 
Comparison oF Proposal To AMERICAN InstiTuTE oF STEEL Construction 
ig. 13 shows the ratio of the cross section areas of columns designed by 
‘the proposed LRFD criteria ia expressed in Eq. 8 (which follows) and the present 
Eq. 9 (also follows). The ratio: is denoted as ‘10 
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assumes | that t the slenderness, I/r of the two o columns’ is 

The comparison is made for different ratios of live and dead load to demonstrate = 
‘the dissimilarity due to varying load factors applied to dead and live loads. 

®. values are presented on the ordinate. 
The | curves were determined by taking the ratio of the _AISC allowable 

‘compressive stress to the stress found by multiplying ®. by F,., in which F,. _ 

is the stress corresponding to R,.. This ratio is corrected by the ratio of factored 
Py in load factor design to the ratio of watactored | loads used in AISC working 


stress design (20). This relationship is developed as shown in | the 


he 
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Fe 
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a 
| The proposed design criteria will provide some economy over the entire range 
of slenderness when the live and dead loads are equal (L/D = 1. 0). ‘When 
live load is greater than dead load (L/D > 1.0), the proposed criteria . may — 
; - Tequire larger columns. If there is only live load, this difference in column — 


= 
= ib | 
ad FIG. 13.—Comparison of Column Areas: Proposed LRFD and AISC 
tm 
in which F,, = R,./A,. Area of column designed by AISC is: 1 peed 
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nlike many structural elements, the strength of steel compression members 
is significantly affected by several factors not quantifiable by the designer. 
2. All column design equations are empirical in that they -Tepresent a range 


of column strengths. However, they cannot be derived from column theory 


as it is presently understood. 


3. Data from 393 tests indicate that there is a large scatter of column strengths “ 
when plotted against slenderness. The variation is due to variables within a 
group of similar columns and to the inherent difference between column types = 
‘and method of manufacture. The variation within one column type appears so 
—_— as large as the variation between different column types examined. 
7 _ 4, Regression analyses indicates that the complex relationships examined did : 
- not account for substantially more of the variation in column strength than > 
_astraight line representation of the dependent variable, nondimensional maximum _ 
ave Tage stress, against the independent variable, -‘nondimensional slenderness. 
_ §, The data base was used to develop regression curves representing column 
_ strength over the entire slenderness range and corresponding ®. values. The 
strength in the elastic range is represented by a hyperbola (Eq. 7). The strength 
‘i the inelastic range is represented by straight lines (Eqs. 5 and es 
6. The proposed design criteria are very similar to the C Canadian Limit State 
criterion. They are slightly less conservative than the criteria proposed by 
Galembos. They provide for columns slightly larger than present AISC if all : 


4 


- designs increasingly less conservative than present AISC until the load is all 
dead load, and the reduction is 15%-20% in column area a aa ia to AISC. Beh 


Aprenoix I. —Recression Equations Using Various Inve INDEPENDENT 


‘The 15 regression equations presented | in the following were developed to — 
ascertain whether the inclusion of more parameters or a parabola, or both, > 
would provide a better fit for column strength than a straight line. The investigation _ 
limited to the inelastic range. 
_Fortyfour column tests were used to develop the following regression equations. —_ 
Only those tests which considered maximum compressive residual stresses (o,), 
shape of the cross section (K), yield stress (F,), and slenderness 
_ The following 15 equations are selected combinations of the foregoing variables. 
The selections are not exhaustive, but are | meant to represent combinations 
leak There are ¢ statistics reported for each coefficient. Generally, 1 t should exceed 
_ 2.0 for the corresponding coefficient to be significantly different from zero. 
7 The coefficient of determination (R’) has been adjusted for degrees of freedom. _ 
7 It is a measure of the portion of the total variability of the data explained 
by the particular equation. The difference between | and R*, (1 — R? ), is 
the unexplained variation. In all 15 cases (1 — mT « ) exceeds 40%. vig ben! wit 
‘This ‘study implies that an increase in complexity of the regression f model on 


pt would not provide any improvement in reducing the unexplained variati 


| 


‘The ¢ statistic for each coefficient given in parentheses and the the first number — 
in each equation equals R? the | lah 


ult 


0.57—* = 0.084 + 0.007 = 


0.44 = =0. 288A + 0.0040, + 1.068 


1.0 ©. 42) 0) 


102 + 0.00020, + 0.134K + 0.012A¢ ,, + 0.837 


Fy 


0.38 ut 

= 0.4400. 0.364n? + 0.00602 + 0.750 8) 
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-=0. 400A 0. 0.427* 


0.52 Ft 0.1370? +0. 
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=0. 00876, 

in which F.,, = maximum average ge column s stress; ‘£F, = laboratory tensile t 
yield strength of column; o,, = maximum residual compressive st stress; and ress; and K 


= sha e constant = 1 if H- ‘sha e and 2 if box- shape. 


Appennix or ResisTANceE FOR Cowumns 


The column resistance factor, ®., is considered by Galambos and Ravindra 
- (10). In essence, it is a measure of reliability within which the true strength 
of a particular element should lie. The mathematics of this determination is 
not pure, so the result should be re-examined against actual test data. In the 

LRFD model, ® is assumed to be a function of the designer’s knowledge about = 

_ the member behavior (professional factor, P), its material properties (material 


factor, M), and the member’s fabrication tolerances and accuracy of manufacture 


_ If the three factors (P, M, and F) are , multiplicative i in the strength relationship, 


“such as R = f(P) x f(M) x f(F), the total coefficient of variation, Vas can 


a 


&§ 


which V = coefficient of variation = = 
by the standard error of the estimate; peo = mean value; P = professional factor 
a error which includes engineering assumptions not considered i in the equation, | 
such as crookedness and residual stress; M = material factor of error which 
includes variation in yield strength and nate modulus; F = fabrication facto factor 


without some corrections. . For example, “yield strength of the shape is very” 

gs in short columns, and of little consequence in slender columns. The 
elastic column strength, in fact, is independent of yield strength. If the error 
due to material properties is a function only of yield 
will be nil for elastic columns. In Eq. 6, , R,. = (1.3 — 0.S7A)A, QF}, 


find that the constant term contains , and the term ‘with lambda 
_ F. Thus the resistance, R,., is related to F) in this range. However, in the 
elastic range, where R,. = )A, OF), ‘cancels, and thus does not affect 


segregation of the three factors in Eq. 26 is not with the 


‘ fabrication such as initial crookedness, straightening methods, and type of <a % 
_ The professional factor also treats variation in material properties such as elastic 

; modulus, which is assumed constant. Yield stress is also included in the variability 
of oo test data ta in that F that F, is used to (0 determine the slenderness, De and th 


q 
= 
| 
4 
| 
pwever the mi) inlic ive accumnti of columns | 


667, 


strength. using ing the measured yield st stress to 

_ the data, much of the effect of yield stress may be removed from the data 

: Set. However, it is known that the nondimensionalization process introduces — 
_ ‘some | error which is biased in an unknown manner. The error may increase 
or reduce the scatter of the nondimensional test data. 
4 _ The addition of a coefficient of variation for fabrication and one for mie 
_ properties may be double counting, in that the most important variables affecting — 
column strength are evidenced in the variety of columns composing the data 


base. ‘the same Wy: and values that w were used by Galambos q 
on 


which is not for in the data ‘because measured colamal dimensions 
“were used to determine area. The V,, term accounts for the variation in yi 
‘stress which has remained unaccounted for previously. > 
Since the equation is fitted to the data, the average value of (P. (P 
is 1.0. The coefficients of variation for m: material and 1 fabrication h agin been n reported 


Ref. 10 as follows: Vy 
Using the weighted COV determined in Table 3, we find Vee 0.10 for the 
= general formulation given’ in Ref. 10 for the resistance factor, ®, is - 


written for or compression members as aT: 


exp (—aB Vz). 
in which a = 0. 55; and 6B = 3.0 (10). First, Vy is found using Eq. 2. For 
the inelastic range, V, = -V.0.10 + 0.10* + 0.05* and V, = 0.15. Then, ® 


‘can be determined using Eq. 27: ®, = 1 x exp (—0.55 x 3.0 x 0.15) = 0.78. = 


_ For the elastic range, Vz = V 0.! 067 + 0.10? + 0.05? = 0.13; ®. = 1x exp 
(-0.55 x 3.0 0.13) = 0.82. 
strength in the elastic range is predicted by Eq. 
A, QF}. If we increase the coefficient such that Eq. 2 becomes R,. 
lhe ‘OF’, , or the Euler curve, the term becomes 0. 937/1. 0. 
“order to maintain the position of the ®.R,. curve in the elastic region of 
®. must be reduced 0. 937/10 > 0 x 0 0. 82 = 0.78. Thus, ‘is 
over the entire range of slenderness. “som 
Assume that the effect of yield strength decreases linearly with pote 
maximum of 7% ath =Otonilath=1.5. 


in which Fym = mean yield stress, in kips per square inch; F’ = minimum 


Specified vila stress, in kips per square i Including the yield si stress. correction 
| and 3 can 


P. 
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nat 
won 1981 Ss 


= 1.28 0.535A — 0.060 + 


0.59% + 0.0257; ; <1.5 


oints vt the relevant range F. /F‘, = 1.28 at A=0 and F,,,/F’, = 0.44 at 


= 1.5 is an appropriate linear replacement. The slope of this line is (1. 28 


— 0.44)/(1.5 — 0.0) = 0.563. The regression line after adjustment for yield 
in the inelastic range and and rounding 1 is ‘al 


— 0.574; A<1.5 Gl) 
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he eine symbols are used | in this paper: 
A = area of column; 
C = Canadian Standards Association (CSA) $16.1 1- 74 column we ah 

E = modulus of elasticity; 

e= maximum eras out- of-straightn 


ettective length factor; 


= column length; 
= load on 
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radius of gyration; 
student’ s Statistic, 
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= load actor; 


= ‘slenderness parameter; 


o in standard deviation estimated by standard error of estimate; and : 


gross; 


ym = yield s stress. ied 
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DIVISION 


1 ‘For Housinc* 


a Egypt has had a long and illustrious history of construction, , stretching back 

at least 5,000 yr to the great era of pyramid building, and encompassing the 

equally great temples and tombs of later periods. That construction was based 

upon local indigenous materials, including stone, soil, silt-based brick, fibers, 
reeds, and gypsum, with its crystalline counterpart, alabaster. 

‘The rapidly- expanding Egyptian population, accompanied by large- -scale shifts, 

f ‘such as rural-urban migration, is placing massive demands upon the —— 

¢ Lt building industry, especially upon its ability to produce housing. This, in turn, 

_ requires immense amounts of materials which must be low in cost, , and therefore, 
they must | be largely, i ‘if not wholly, indigenous. Recognizing these needs, the 
= Technology ‘Adaptation Program carried out by the joint collaboration of alll 

University and the Massachusetts Institute of Technology under the auspices . 
_ of the Agency for International Development of the United States Department er 
of State has concentrated on indig2nous materials in the portion of its program, 
_ devoted to housing to meet the needs of the Egyptian Ministry of Housing 
the Ministry of Development and ‘New Communities. tga) 


- Among the materials considered are soil and soil stabilization, low-cost cement 
lightweight aggregates, and concrete, brick, fibers, and gypsum. Gypsum is a 
an attractive, now little-used indigenous material. This paper examines strength § 7 
and related properties of gypsum | blocks « and nd panels with various additives a 7 
on chosen because, after the loss of the Sinai Peninsula its 
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7 only, and might therefore help to relieve the pressure on tight supplies of portland 
Cement concrete products and dwindling supplies of the traditional red ‘brick | 


paper presents some aspects of the work o on n gypsum 2 as an item of possible 
value to others interested in materials, especially indigenous materials, for 
_ construction in general, and housing in particular, = ea. 


An extensive experimental program was undertaken at | Cairo” University to 


study the different characteristics and potentialities of locally-produced gypsum 
_ for use in various components for walls and roofs in housing units of different 
TABLE 1.—Physical and Mechanical Characteristics of Locally-Produced Gypsum" 


 -‘Testresults 
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with 0.1% (ASTM-CAT: C472) | Water / plaster ratio by 
Setting time (water /gypsum = 
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Ot 
48, no o keratin) nen 42.21 kg/cm? ad) 


= All results are based on the average of five tests; (maximum deviation within +10% 
into three main p parts dealing w with p plain gypsum, 
gypsum masonry blocks, and reinforced gypsum components. This paper presents" 
the work on plain gypsum and gypsum blocks. 
Program on Plain Gypsum.—This part of the program was undertaken on — 
~ locally- -produced plain gypsum to determine its physical and mechanical proper- 
i ties, and assess its resulting modified characteristics as affected by the percentage 
(of mixing water, the degree of compaction, the method of drying, the amount 
of setting retarder, the addition of sand, the presence of sawdust, and the 
amount of sodium bicarbonate added to foam gypsum. 
Physical and Mechanical Characteristics. —Locally- produced hemihydrate beta 
= was used in all tests, which } complied | with standards of the American 
Society for Testing and Materials. This included: fineness, time of setting; and 


A 


of high-grade gypsum. It appeared possible, therefore, that gypsum might well | 

: he nsed for more diverse huildine nurnoses than its traditional role as niaster _ 7 
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FIG. 3. _—Effect of Sand /Gypsum Ratio on: (a) Compressive | Strength of Gypsum; 
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were left, after casting, in moist air for 24 h, then placed in an over at 104° 4 
‘(40° C) and tested after seven days, when the weight had become constant. 


Percentage of Mixing Water. _—The effects of varying the water to “gypsum _ 
‘Tatio, degree of compaction, m, and method of ‘drying on the resulting s¢ seven-day » 
_ compressive strength of gypsum without keratin retarder were studied for three a 
groups: A (not vibrated), B (vibrated), and C (vibrated and oven-dried), as 
shown in in snes i. The variations in density and relative workability of vibrated 
Sawdust Passing Sieve ‘Sawdust Passing Sieve inch, 
W/G Variable, K/G = 0.1% 7 W/G Variable, Keratin = 0. 1% 
wi 
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ajo 
PERCENTAGE OF soDIUM BICARBONATE - OF sop1uM BICARBONATE - 
(& by weight of gypsum) (% by weight of gypsum) 


of Foamed Gypsum; (b) Density of Foamed Gypsum —— 


specimens as affected by the water to gypsum wn ratio are in 


and l(c), respectively. Workability Curing ‘mixing was judged by the relative AR 
ease of mixing and of filling the molds. 
, g@ Amount of Setting Retarder (Keratin).—Keratin is commonly y employed to 
retard the setting of gypsum. Different mixes were prepared with percentages _ = 
of keratin ranging between 0.05% and 0.20% by weight of gypsum. The resulting : 
setting times and compressive Strengths ¢ after seven days are plotted in Figs. =~ 
2(a) and 2(b), respectively, —2 08 
Effect of Adding Sand. d.—Sand passing a No. 14 sieve was added to gypsum 
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«FIG. 7.—Cross Section of 29 cm x 29 cm x 10 cm Cored Wall Block '§ 


TABLE 2. —Mix of Solid Blocks 


‘Filler material 
or additive 


Sodium bicar- 
bonate and 
‘sawdust 
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or 


‘Sawdust 


sodium bi- 
carbonate 


“Sawdust to gypsum ratio is by volume. 
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Block number | ore 
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Blocknumber | G | W | S | ate, N 
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TABLE 3.—Mix Proportions of Tested Cored Gypsum Blocks 
1 | | — {82 
1 | 00 | — 0.006 
1 jos | — [a ooo 
“Sawdust to gypsum ratiois by volume. 


without | keratin in ratios by “weight aun tan 0.5-1.5. The 
compressive strengths after seven days and the densities are plotted in Figs. 

4 3(a) and 3(b), respectively. bar. 
is Effect of Adding Sawdust.—Sawdust passing a 3/8-in. (0.95-cm) sieve was — 
4 added to gypsum with 0.1% keratin ratios by volume ranging from 0.5-2.0. 

7 The resulting densities and compressive strengths after seven days ¢ are plotted = 

in Figs. 4a) and 4(b), respectively. Varying wa My to gypsum ratios were adopted _ 

to attain excellent workability for all mixes. 

Foamed Gypsum.—Gypsum might be to reduce weight and increase 
thermal insulation. For foaming gypsum, sodium bicarbonate was added to the | 

dry gypsum in percentages by weight ranging between 0.2-5.0. Keratin was a 
added to the dry mixture as 0.1% by weight of gypsum, and the water oo 
gypsum ratio by weight was 0.45 for all mixes. Compression test cubes were — 
cast without mechanical vibration to avoid losing a part of the foam. 088 | oT 

_ The cast cubes were tested after 14 days in laboratory air when their <= 

become constant. The effects of varying the percentage of sodium Cesbensts 


_ Compressive St rength is 
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nd 
FIG. 8.—Effect of Sand/Gypsum Ratio on: (a) Compressive Strength of Solid and 
- Cored Gypsum Blocks; - Dry Weight of Solid and Cored Gypsum Blocks 


| 
Gypsum Masonry Blocks.—This part of the program was intended to develop 
_ masonry wall units made of unreinforce¢ sh a filler might be : 
“Tl 


added, if needed, to decrease weight, cost, or both. 
i The dimensions of wall units had to | give a a thickness allowing the grip of 
: 2 “the hand, and at the same time give a reasonable stiffness and impact resistance 
f = to the plain gypsum wall unit. The overall size had to be as large as possible 
to decrease handling labor during construction and the quantity of mortar required. : 
Yet the size was limited by weight, allowing the unit to be lifted by one or 
two hands. Taking the preceding into account, a wall unit with overall all dimensions, 
SS 42 in. x 11.42 in. x 3.94 in. (29 cm. x 29cm x 10 cm) was chosen, satisfying © 
- module of 11.81 in. (30 cm) and allowing for 0.39-in. (l-cm) thick mortar 
_ Manufacturing of Gypsum Wall Units.—Gypsum wall units were case in 
aluminum forms consisting of a base and four sides connected by quick-release 
bolts. Three types of forms were , used. The first gave a tongue-and- groove 
type block as shown in Fig. 6; the second gave a block with flat sides; and 
the third had flat sides containing circular openings in two opposite sides through ~ 
which plastic pipes could | be inserted and removed, giving a cored gypsum — 


me Sand, | sawdust, and ‘sodium bicarbonate | were tried a as fillers or additives t to 


 (eoratia) to gypsum was 0.1% by weight. Sand passed a No. 
passed a 3/8-in. (0.95-cm) sieve. 
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f 
‘The cast gypsum blocks could be demolded 30 min-45 min after casting. 
ul ‘Test Procedure.—The gypsum blocks were weighed just afer demolding and 
placed in ‘Telatively high h humidity for 24 h, then removed to the open air and © 
- weighed every day until completely dry. All types of solid and cored gypsum 
3 blocks were then tested in compression up to ultimate failure. Blocks with 
_ tongue and groove were first prepared by trimming the protruding tongue, and 
filling the groove on the opposite side with appropriate mortar. 
2 ‘a Test Results.—The effects of type and percentage of additive on the resulting — 
ESE strength and dry weight of solid and cored ‘gypsum blocks are 
ny in Figs. 8-10. The resulting material cost of solid and cored eypsum 
block affected the added of sand or sawdust is plotted in 


PERCENTAGE OF SODIUM BICARBONATE by of gy psum) 
FIG. 10. —Effect of P Percentage of Sodium Bicarbonate on: (a) Compressive Strength 
of Foamed Gypsum Solid Block; (b) Dry Weight of Foamed Gypsum Solid Blocks. 


Figs. 11 and 12. The effect of the existence of cores on the drying time needed 
for pure gypsum blocks and blocks of gypsum mixed with s sand is shown in 
‘The tests indicate that cored blocks dry eight days-ten | ten days earlier than’ ; 
equivalent solid blocks, and that higher percentages of filler or additive, especially 
7 Recommendations. —Based on these results, the following types of blocks 
 -* Solid blocks with a sand filler, where the sand to gypsum ratio by aime  &§ 
is 2.8, are recommended for load- -bearing walls; as their 


' 
4 PERCENTAGE OF SODIUM BICARBONATE (% by weight of gypsum) 
| 


is about 639 psi (45 kg/ cm. and their weight is about 32 Ib/ block | 
(14.50 kg/block), whereas the cost is about 0.063 L.E./block ($0.09/block). © 
4 The required water to gypsum ratio is 0.77 by weight to give excellent workability. ' 
_ The value of the Egyptian pound was about $1.45 at the time this report was _ 
| Cored blocks with sawdust as a filler and sawdust | to gy gypsum ratio of 
0.8 by volume may be used as nonbearing wall units as their compressive strength  « 
is only about 213 psi (15 kg/cm’; 1,469 kPa). The corresponding weight is a ot 
about 15.4 Ib/block (7.0 kg/block), and the cost is nearly 0.078 L.E./block 
4 ($0.11/block). The required water to gypsum ratio by weight is 0.66 to give 
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FIG. 11.—Effect of Sand/Gypsum Ratio on Cost of Solid and Cored Gypsum Bi Blocks 
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3. Whenever weight is not a problem, a lower cost can be obtained by ving 
ae blocks for nonbearing walls with a sand to gypsum ratio ) by weight s 
2 2.5 giving an approximate cost per block of only 0.05 L. E. ($0.07). The 
corresponding weight is 24.2 Ib/block (11.0 kg/block), and the compressive 
strength about 213 psi (15 kg/cm’; 1,469 kPa). The water to gypsum ratio 
¥ The cost of all the preceding recommended blocks compares quite well with 
the cost of an equivalent wall area of traditional red bricks, which would | be 
not less than 0.15 L.E. ($0.21). 


— for Suggested Gypsum Blocks.—In an effort to ‘study the type of - 
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_ as 0.1% by weight of gypsum. The resulting compressive and tensile strengths _ 
3 after seven days were recorded for the different mortar mixes. The water to 
gypsum ratio by weight was varied for different mixes to attain excellent — 

_ workability. The compressive strengths for the different mortar mixes are plotted 

in Fig. 15. Test results indicate that a mix Proportion of gypsum:sand: ‘lime — 

bya volume of 1:1:0. .2 and a ratio | of w water to gypsum by a , weight of 0. 70 

~ would give ac compressive strength after seven days of 1,107 psi (78 kg/cm* : 

7,634 kPa), and a corresponding tensile strength of about 255 psi (ig kg/cm 5° 

‘which is is quite for use as a mortar. 
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- recommended type of mortar, tested on six bond samples ; as shown in Fig. a 

16, was found to be 64 psi (4.5 kg/cm?; 441 kPa), which is satisfactory. 

_ Group Action of Suggested Gypsum Blocks.—In an effort to study ne 
action of the developed gypsum blocks, a series of test walls were built using 
& suggested blocks and recommended mortar, and tested in compression and — 
flexure. Six test walls were built using load-bearing gypsum blocks with plane ~ 
SS to be tested in pure compression, one in pure flexure, and 
the rest in combined compression and flexure with ratios of applied compression 
b ultimate pure compression, (P/P,), of 0.71, 0.41, and 0.28. Three test walls 
_ were also built using load-bearing gypsum block, with tongue-and-groove connec- — 
tions, with one to be tested in pure compression, one in pure flexure, and a 
the third in } combined compression and flexure with (P/P, ) of 0.50. A test 


was s tested i in n pure flexure. All of the ten aforementioned test walls w were 3. 94 
in. (10 cm) thick, 23.62 in. (60 cm) wide, and 47.29 in. (120 cm) high. To 


_ a mortar suitable for the suggested gypsum blocks, and to avoid the use of cement © P 
martare tn cave cement far hetter unease caverel auncum mives ware 
| 
FIG. 16.—Bond S h Test Samples 
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4. —TostF Results of Loaded 
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pe 
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blocks (groove and 


Cored gypsum (1 (nonbear- 


Traditional, red bricks 


bricks 


b "provide: for - comparison, a test wall of sand-lime bricks and two test walls of 
traditional red bricks were built with the same dimensions as the gypsum test | 
_ wails, except for the thickness which was 4. 4.72, in. (12 cm), and pasate 
oe full description of the different test walls and their corresponding recorded 
-- results is given in Table 4, while the characteristics of the types of blocks sd 
in building the test walls are givenin Table5. on 
_ A test setup shown in Fig. 17 was specially constructed in this investigation * 
to ber mounted on a Amsler Compression Machine to givea 
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_ EGYPTIAN GYPSUM 


in Combined Compression and Flexure ie > 
Ultimate | Ultimate | 
Maximum bending "stress, (in stress, tin 
compressive moment, | kilograms 4 kilograms 
load, P (in “in kilogram per square per square 
| 


(11) 
Vertical crack through 


Vertical crack in upper a 
row only 
Horizontal failure a 
midheight 
Horizontal failure at 
_midheight 
| Horizontal failure 
midheight 
Horizontal split ; 
block and mortar 
| Vertical crack in — per 
rowonly 
Horizontal at 
midheight 


midheight between 
block and mortar — 
Horizontal split at 
midheight bet ween 
block and mortar 
Vertical il crack through 
wall 
2.89 Horizontal split at 
block and mortar 


For walls tested under combined compression and flexure, the total predeter- : 
mined compression load, (P), was applied first; then the flexural load was applied — 
in increments up to ultimate failure. te 
interaction diagram between flexural and ‘compressive stresses was expel experi- 
= mentally developed for load-bearing gypsum blocks with plane connections, - 
s is shown in Fig. 18. Such an interaction diagram gives the nore 


combinations of Senueal and 1 compressive stresses causing failure of the loaded 


« 


a 


program ‘on test walls indicated pure compressive ‘strength of 
load- bearing gypsum walls is nearly twice that of the traditional red brick walls, 


and is is is greater by at least 30% | than that of sand-lime brick walls. — 


TABLE 5. —Characteristics of Block Types Used in Tests* 


volume,| Ultimate 
"Average Overall Dimen- 


in kilo- sions, in centimeters 


per square 
Type of block grams" Length Wi 
Cored nonbearing gyp- 
sum block 
Traditional red brick 
Sand-lime brick 


_ Members of Supporting Frame 


"govtes to Supporting Frame 


FIG. 17. BaP Moy Sketch of Loading Test aati and Supporting Frame for Flexural = 
_ The pure flexural strengths of the different bearing and nonbearing gypsum = 
walls were found to be nearly equal and greater by about 60% than the pure i 
"flexural — of a corresponding traditional red brick wall. On the other 
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EGYPTIAN GYPSUM 687 
pry gypsum walls with both plane and tongue- -and-groove connections seemed 
A to have no particular difference in behavior under load. Tongue-and-groove 


Bearing Blocks with Plane ol of | ve 


_ blocks, however, required more care during construction and consumed. a Hlonger 
_ time and a a larger amount of mortar than the blocks with plane flat sides. SP ions a ot 


Based on the tos tests, the recommendations can be made. 


tas 


tountar while the comk vaced (1). 
Ys Solid gypsum sand | blocks should be carefully considered for load-bearing 


walls and partitions carrying moderate loads at least equal to those tradition nally 
2g _ supported by conventional sand-lime and red bricks. = ees . Few have 
~ 2. Cored gypsum blocks with sawdust filler where light weight is important, — 
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NAL OF TH E 
. STRUCTU RAL DIVISION. 


By Thomas Brown, ' A. 

_ With the current emphasis on structural optimization and the trend to limit 
states design, it is essential that we obtain a proper understanding of the instability 

7 “(and the factors influencing it) of ‘nonprismatic structural members. While a 
2 number of | both theoretical and experimental investigations have been carried 

— out, it is apparent that the buckling behavior of nonprismatic structural members 


under all possible loads and load combinations is not yet fully documented 


The behavior of tapered columns acted on by axial load or combined axial 
_ load and end moments has been the "subject of a number of investigations © 
(4,6,7,10, 13,16). These have resulted in interaction curves similar to the one 
a used for the design of prismatic members (17). The lateral-torsional instability : 
‘ of doubly-symmetric, structural sections under the action of lateral load or 
end moments, or both, is a complex problem without the added inconvenience © 
- of a variable cross section (1,18). Timoshenko (18) first identified the problem — 
and presented solutions to a number of cases of prismatic beams under various © 
transverse loads. In addition, Timoshenko presented critical loads for narrow 
rectangular cantilevers in which the depth varied (18). Salvadori (16) used the 
: Rayleigh Ritz method to evaluate critical end moments for a variety of moment 
gradients. A number of investigators (2,8,14) have used various methods to 
obtain critical loads or moments _ for a number of particular cases. Few have 
considered the point of application of the load other than at the centroid. Butler — 
4, 5) has presented the results of a series of experimental programs none of 
which consider the case of pure lateral-torsional buckling under transverse load. "i 
7 More recently, Kitipornchai (12) has presented solutions for the elastic stability _ 
simply-supported tapered I-beams including comparisons with 
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results. In this paper, critical loads for beams, both -eupported 
and cantilever, are determined using finite differences. Particular attention is 


paid to the location (with respect to the centroid) of the transverse load. a) 


The characteristic eigenvalue equation for any structural stability problem 4 
may be developed from: (1) Equilibrium considerations; (2) imperfection e equilib- 
rium; and (3) energy considerations or dynamics. ae 

The pre pre-buckling and a -buckling configurations of a ‘donbly-cymaetric 
I-section are shown in Fig. 1. The displacements of the centroid are denoted 


as u. and v.; the dieplacemeats of all all other points on the section are given 


in which the prime denotes differentiation with respect to the longitudinal axis, 


parse: 


; and >, = the warping function. Using the foregoing displacement functions, 


+B (v’, + xB’ 


u 
q 
a 

om 


components: (1) The to the strains acting on the stress 


field immediately prior to buckling; and (2) the strain energy ines by th tee . 


j may. 
oy€ydv + 


n 
For “open thin-walled sections, the following integrals result as 


+ u’)? + El,u”? + EY —2ER 
4 
+ + GAB + dx 


in poten entia ener 
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becomes, in i nota Eqs. 5 and 6, and introducing the work term, W: © 
ow 


ll 5 
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(EI, u"? — 2ER,B” + GKB’? + ET dz + boundary terms (8) 

the for potential ‘energy, rgy, the terms associated with 

Baler buckling about the major (x) axis have been omitted and fourth order 

_ Eq. 8 for the change in potential energy is developed with respect ‘to the 

centroid of the cross section and this accounts for the minor dif! ferences between 
og this equation and the corresponding equation in Ref. 1 and Eq. 10 | of — 

Considering the variation ad 8, we 2 


aed 


ae 


[- 
+——-B’ | dz...3; 2! =0; 2 )dv 
=F total change 
| 


APRIL 1981 


(EI, -E ER,p” y" (MB)” +(P, u Ae =0. 


‘The result, then, is two simultaneous fourth onder differential equations. Unter : 
Certain conditions, , these equations may be ‘uncoupled; for 


— 


result in 


wae the equations become 


ay 


boundary conditions, leads to the equation 


: Depending on the beam type, M will be so some function of z and the load. 


can be approximated using central finite (9), leading to” 


Eq. 12 be twice and, with the of 
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TORSIONAL BUCKLING 


This equ: iia is of pon form and therefore. does no not lend itself to rnd 4 
~ eigenvalue extraction techniques. It may, however, be solved directly if we 
es Premultiplying Eq. 17 by [4|™ ana combining Eqs. 17 and 18 results in 


from which ch eigenvalues extracted using any conventional eigenvalue 


program. 


RITICAL Bucxunc Loans 


0. 727 0.783 
0.835 
0.849 0.882 
0.903 | 23.37 | 0.924 
0.952 0.963, 
16.44(5,) | 1.0 | 25.28(,) | 1.0 | 


width and thickness and the web thickness rer remaining constant. The: section 
used for all examples was an idealized I-section with a maximum depth of = 
24 in. (610 mm). Flange width and thickness were 6 in. (152 mm) and 0.5_ 


(13 mm) respectively, and the thickness was 0.375 in. (9.5 mm). 
Critical loads are given in form = os teh Che 


L’ 


J the case of nonprismatic simply-supported beams 
. __ and nonprismatic cantilevers. In all cases only the beam depth varies; the flange 
ABLE 7, "or mene unter End 
TABLE 1.—Coefficients for Simply-Supported Beams under Central Concentrated 
Taper ratio a 
— 
0.167 
9667 
: 
7 
| 
| 
| 
rondiuions 


expected, the top results i in significantly reduced cri ical loads 
as compared with the other load positions. There is not, however, the large _ = 
reduction in critical load with the rather extreme case of an end depth of 4 

Table 2 summarizes the coefficients y, for critical distributed loads, 
2 

acting on 20-ft (6.10-m) simply supported beams. 


ration 

0.834 
067 | 2.34 7 


= 


ei , 7 


¢ 


a, 
| 
in Fig. 3, the boundary conditions used 
| 
= 


i 
- Table 3 indicates that a 2 thin | in the free end depth of a cantilever leads 


to an increase in the critical end load [Fig. 3(a)]. This is due to the very 
significant effect of the top-flange load. This is demonstrated by Table 4 which 
shows the effect of varying the point of application of the load from the centroid 
(a/d = = to the top (a/d = 0. 5) of a prismatic canti cantilews er [Fig. 


TABLE 3. —Coefficient Y, for Cantilever Beams under End Concentrated Loads 


F,) 


‘e 


4.—Coefficient Y, for Prismatic Cantilevers under End Concentrated Leeds 
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the corresponding | load at tk at the centroid. 


Central finite differences are used to determine critical ‘buckling loads 


for 


mt included. This effect is shown to be significant in all types of beams but. 
leads to an ene, -, wm decreasing free end depth, of the critical load in 
are due to Omar Al-Amoudi and Abdul-Rahman Khathlan who « 
gut all of the computational work and to John F. Doyle for advice on the 


mathematics. -t 
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= vertical distance between point of application of transverse _ 
depth of I-section; 
modulus of elasticity; ‘a 
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shearmodulusy; 
= principle moments of i inertia; BRN AL : THE 
moment of inertia at deepest section 
K = St. -Verant torsion constant; wall SION: 
k= Sor _ 


= total bending at any point in 


a M, = moment due to eccentricity of axial load; 
concentrated transverse load; 


axialload; 


characteristic 


distributed transverse load; 


R,,R, = SSyo.dA, §§xo.dA; 
w= displacements of any point in in beam; 
We = displacements of centroid; = 
total potential energy during 


work done by external loads during buckling; 


coordinates; 


ratio; 
angle of twist; 
amended ‘eigenvector (Eq. 18); 


= 
radius of gyration; 


additional stress tensor caused by aed 7 
stress tensor prior to and 
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Notes. 
[AUSUGL 
To provide a place within ASCE for publication of ideas not d 


as yet, to the point where they warrant publication as a Proceedings paper in a Journal, the | 
7 publication of Technical Notes was ee 7 the Board of ines on October 16-18, 


1967, under the following guidelines: _ i 


1. An original manuscript and two copies a are to be submitted to the steaiiiee of — 
; and Professional Publications, ASCE, 345 East 47th Street, New York, UN. -Y., 10017, along = 
with a request by the author that it be considered as a Technical Note. 
» 2. The two copies will be sent to an appropriate Technical Division or Council for review. | 
3. If the Division or Council approves the contiution for publication, it shall be returned ; 
4. The technical publications staff will prepare the material for use in the earliest possible 
issue of the Journal, after proper coordination with the author. __ er Ore a 
. Each Technical Note is ‘not to exceed 4 pages in the Journal. As an approximation, 


. The Technical Notes will be grouped in a special section of each Journal. 
: fone retrieval abstracts and key words will be unnecessary for Technical Notes. a 
. The final date on which a Discussion should reach the aes is _ as a footnote — 
‘each Technical Note. 


5 1. Titles must have a length not exceeding 50 characters and spaces. == oo 


he 2. The author’s full name, Society membership grade, and a i footnote reference stating present 


3. The manuscript is to be submitted as an copy two duplicates) that is typed 

_ double- -spaced on one side of 8-1 /2-in. (220-mm) by 11-in. (280-mm) white bond paper. a 
4 All mathematics must be typewritten and special symbols must be properly identified. 


The letter symbols used must be defined where they first appear, in figures or text, and arranged 
—& Standard definitions and symbols must be used. Reference must be made to the lists 


published by the American National Standards Institute and to the Authors’ Guide to the 


Tables must be double- spaced (an ribbon copy and two duplicate copies) 


be submitted as glossy prints. Explanations and descriptions must be made within the text 
8. References cited in text must be typed at t the end of the Technical Note in . alphabetical 
9, Dual units, i.e., U.S. Customary fe followed bys SI (International System) units in n parentheses, 


"should be used throughout the he paper. a RA 
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IN HIGHER Moves 


B Ahsan Kareem, M. ASCE 

tall buildings i is found to y 
in the fundamental sway modes. The higher modes of vibration have an 

insignificant contribution to the displacement response; however their influence — 
on the higher derivatives of response, i.e., acceleration and jerk, is quite significant 
® ,4). The human biodynamic response to motion is ‘generally influenced by 
the higher derivatives of displacement (4). Therefore, occupant comfort in tall 
buildings can be affected by underestimating the dynamic response as a conse- 

quence of neglecting the higher modes of vibration, 
_ Simiu (5) has shown that the contribution of the higher modes to the ro aa 
-mean square (RMS) alongwind acceleration is on the order of 5%-20%. - Fujimoto, z 
Ohkuma, and Amano (1) utilized a four-mass aeroelastic model to. study the 
= response of a prototype building. They : similarly concluded that the 
contribution of the higher modes to the RMS value of the acceleration response — 
a a was 10% of the response in the fundamental mode; waenens the displacement 

< response was not influenced by the higher modes. = corte) 
For high wind speeds, the crosswind response of a uniform cross-sectional 

tall building is generally higher than the alongwind response (2,3). In this note, J 

a simple method of computing the contribution of higher modes to the crosswind 


_ Tesponse is presented, and the veuulis are compared with those of a detailed 
| 


= It can be illustrated from a dimensional argument ‘that the mean square 
acceleration for a lightly-damped structure with narrow-band response venensilll 
teristics is proportional to f, S(/,)/é,, in which f, = the mth natural frequency; 
€, = the structural damping ratio in the nth mode; and SU.) = the ‘power 
spectrum of the wind- -loading function evaluated at the nth frequency. 
_ The power spectrum of wind loads can be expressed as being proportional as being proportional — 
"Asst. Prof. of Civ. Engrg., Cullen Coll. of Engrg., Univ. of Houston, Central Campus, Houston, to oy 
Note.—Discussion open until September 1, 1981. To extend the closing one 
a written request must be filed with the Manager of Technical and Professional Publications, 
_ ASCE. Manuscript was submitted for review for possible publication on October 8, 1980. 
This paper is part of the Journal of the Structural Division, Proceedings of the American 
} Society of Civil Engineers, ©ASCE, Vol. 107, No. ST4, April, 1981. ISSN 0044-_ 
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6in which 
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the building geometry (2). For the crosswind forcing function (3), the value 
of a is to 5. scceleration Fesponse is given 


(1) 

between the bp and the first — contribution to the acceleration re 


In the higher the relatively high flexural shear of 
the structural system result in higher damping than occurs in the fundamental 
mode. Yokoo and (6) established an approximate expsession for 


‘in the fundamental mode in terms s of the a: associated frequencies 


the R RMS value of the acceleration The ratio 


in whieh C= = approximately equal to 0.38. Eq. 3 into 


"contributions to various components of response is given by hie 


f in which r represents various derivatives of the displacement response; ¢.g., 


- (the third derivative of displacement) is represented by r equal to three. 
_ The wind-excited response of a tall building modeled as a multidegree-of- 


freedom system can be estimated using methods of random vibration theory — 


(2). The mean a square acceleration at the ith level is given =. = 


(G,(f/)] = (14, [6] (i2m f)}}, 


| 
|. 
| 
| : | 
| 


= total number of modes. tar 
It will be assumed that the cross terms in Eq. 7 are relatively small and 
that the forcing function varies slowly in the proximity of the building frequencies. 
The forcing function thus can be idealized by a white noise process. This reduces y 
Eq. 6 to a simplified and approximate relationship) 


and the. mass in the » nth mode, respectively. A m more general expression for 
mean value of various components of response is given by 


in which ), m the generalized power spectrum of the function, 


a tne 100 ft Gl m) square in plan, and 600 ft (183 m) ) tall was as modeled 


a five-degree-of-freedom system (2). The building characteristics are given 
in : Table 1, and the kinematic model is shown in Fig. 1. Stemi | 


contribution of the crosswind response in the higher modes was computed 


‘TABLE 1.—Building Characteristics” 


—0.6 —0.8 


| 


— sm TECHNICAL NOTES 
| 
. 
, 
| 

“Building density = 12 | 


Ea. at for of The results are listed in 
Table 2. The displacement, acceleration, and jerk components of response were 
also computed using Eq. 9. The results are listed in Table 3, in which the 
numbers in the parentheses represent the contribution of each mode as a 
_ percentage of the first mode. The building response in the second mode is” 
"significantly affected by the acceleration and jerk components. _ This trend, 


_ however, tends to diminish in the succeeding higher modes. The estimates 


4 

4 

«4 


computed from Eq. 5 are Po ee then the vated’ obtained from Eq. 9 for the 
second mode but tend to be higher for the succeeding modes (see 1 Tables 2 
and 3). This trend can be attributed to the changes in the value of « for the | 
crosswind forcing function with an increase in the reduced frequency > Se 
_ For the example in this examination, the value of a near the first two frequencies 
is lower than 5. 3; however for the third and higher frequencies, a tends to 


-¥ 


| 
| 


exceed 5.5 value used in the present analysis. Therefore idiaaaaaiiines 
appropriate segmental values of a, one could achieve good « estimates of ‘Fesponse 
The u use 2 of the ‘simple e expression gives by Eq. 5 for f preliminary | design purposes dq 
would allow early of structural requirements needed to resist 


TABLE | 2.— -—Contribution o of t Higher Modes /First | Mode, ase percentage 


=. 


TABLE 3.- Modes Contribution, for Level 


0.0002 
as) ax _@. 


=... 


0.0272 0.0058 


“Numbers i represent of first mode. 


the oscillatory building response, and appraisal of the need for a more detailed 


‘The mates herein | suggest and "stability criteria 

- for the design of a tall building can be sufficiently satisfied by a fundamental 
= mode analysis since displacement response is completely dominated by the first 
mode. However, the assessment of serviceability and occupant comfort criteria 
requires | ‘the: inclusion of higher modes of vibration in the dynamic analysis — 


& 
0.219 0.0000 (0.000014 
0317 0.0412 | 0.0054 0.0078 
ner ie oom 
RMS Jerk, in feet per cuble second = 
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-....- or on any paper presented at any Specialty Conference or other meeting, the Proceedings 
: f which have been published by ASCE. Discussion of a paper /technical note is open to - 
any 
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Closure by Allan P. Mann‘ and Linden J. Morris* 


A difficulty that arises in comparing wing what at first appears to be similar research 
evidence from diffi erent sources is that there might i in fact be significant differences 
in the purposes of the separate research projects. This in itself may be of | 
i little consequence provided other researchers and practicing engineers are made 


fully aware of the aims and limitations of the particular research projects being _ 


_ reviewed. For then the potential users of the research information can realistically 
evaluate what may be different recommendations for the same design problem. 
It is, therefore, important that the writers reiterate the aim of their own research ~ - 
_ and recommendations—which is to provide a design method for extended 
end- plate connections by which the beam moment reaches its plastic moment | 
~~ while the connection sustains only limited distortions which neither 
inhibit the joint strength nor its rotation capacity. The recommendations are 


in the context of having four bolts grouped around the beam tension flange. 


_ Connections proportioned by the writers’ recommendations will have the 


capability of providing adequate rotation capacity while still sustaining the plastic ~ 


moment capacity of the beam. This will allow the redistribution of moments ~ 
in the frame as required for plastic analysis. On the other hand it appears, 
though not specifically stated in so many words, that the discusser b bases his 
proposals < on the performance of end- plates subject only to service loading, 
with an implied reserve of strength, i.e. 1. 7x maximum elastic —— 
| 7 So working load. His test results (28) indicate that in most cases the plastic 
capacity of the beam was just reached. However, the M/M Tatios achieved 
represent a necessary but insufficient requirement for connection design. -Ade- 
quate rotation capacity is just as important, but very little evidence of 
4 moment-rotation characteristics for the discusser’s tests is given. Furthermore, — 
_ itis stated (28) that the end-plate thicknesses produced by the discusser’s method 
- would frequently be accompanied by some reduction i in connection rigidity = 
under circumstances where deformations are critical, the influence of such | 
increased flexibility may need to be reviewed. It should also be remembered — 


flange is stiff or adequately stiffened, though neither term is defined, Pe ; 
the writers’ procedure covers other practical cases where the flange thickness — 


is less than or approximately equal to the end-plate thickness. The research | 


1979, by Allan P. Mann and Linden J. Morris Paper 


*Lect., Dept. of Engrg., Univ. of Manchester, Labs., Oxford Rd., 


| 
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by Packer and Morris (14) showed that it was Perr deformation -—_— 
; the column flange and end-plate that governed whether or not the moment-rotation - 
_ characteristic of a connection was satisfactory. Fig. 13 shows the effect of 

_ different column flange thicknesses on a connection’s moment-rotation behavior. ‘y 
If f only the attainment of strength was considered then it could | be said that _ 
specimens J1, J2, and J3 were all satisfactory having achieved moment ratios, 5 


_ M/M,, of 1.21, 1.09, and 0.98; however when rotation is also a — 


Retetion inradions x 10 
FIG. 13 ciliated Rotation Characteristics of Connections /1 to J3 


att 


‘ 


then en only J behavior implicit in the: writers’ method. 
Similar curves can be produced for decreasing end-plate thicknesses. Indeed, | 
the discusser refers to tests by Beall (22) which seem to confirm that the connection 

_ moment capacity decreased with a reduction in end-plate thickness. By reducing 
the column flange or end-plate thickness below that recommended, the effective 
““moment”’ ’ capacity of the bolt group, and thus that of the connection, Pomscentin 


fora given rotation requirement. mets 


a! 
> 
— | 
| 
= 
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a not be considered in design. Certainly the writers would agree with that opined 
iy if only allowable stress conditions were being considered. Figs. 14 and 15 support 
this argument as any prying action is nominal at working load level. But the 
oe writers would contend that the behavior of the bolts between service and limiting 
bad plastic loads has an important bearing on the design of the bolts even in allowable 
_ stress design where there is an implied reserve of strength of 1.7 x working — 
load conditions. The diagrams show that there is a definite enhancement in 
the predicted direct bolt load (F,/4) which if not allowed for in design ll 
t result in bolt failure. As Agerskov noted in a previous discussion (32) in the _ 
case of specimen J5 ‘Fig. 14) the apparent ‘Prying force i is some 10% i in excess > 


connection ‘does ‘not in “fact comply with the writers’ recommendations, i 
column flange requires stiffening. Specimen J4 (Fig. 15) represents a stiffened 


Lt... ‘ded 


whe 
16. Tension Bolt Load Behavior for 


_ JS and accordingly the measured bolt forces show that the prying force had 


been reduced to about 33% the estimate of the prying force assumed in the 
design procedure. The value of 33% was chosen as representing a —— 
upper limit on the increased bolt force due to prying action for connection 
designed by the writers’ method. This results in a simple design rule for —. 
_ The discusser agrees that the inner row of bolts will attract more of the 
_ direct load and concludes that these bolts could become somewhat overstressed 
7 at service loading. The discusser then suggests that the possible development 
_ of prying action (!) would be resisted wholly by the understressed outer row 
- of bolts. This does not negate the problem of the overstressed inner bolts! — 
Indeed, these inner bolts in a splice connection have been known to fail = 


flange. As ‘explained in the original paper, in the context of beam- “ll 
connections, the lower bolts carry a larger share of the direct force (F,) but 
as prying action increases the upper bolts attract more of the prying force 


until at limiting plastic load level the total force (direct plus prying) in both | 


7 
— a 
4 
— 
1 
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os inner and outer rows tend to become eqeal: ‘This is by. ‘the 
measured bolt forces indicated in Figs. 14 and 15 for the specimens J5 and 
-? Clearly, if the prying forces are taken into account in proportioning the 
bolts then there should be no problem with overstressing the bolts. wily gra a 
ty The problem of not having any bolt failures is that it is difficult to assess 
the reserve of strength left in the bolts, i.e., is it 1% or 101%? It is always a 
helpful to the researcher to have the occasional bolt failure in order to judge = = 
the performance of the bolts. The discusser suggests that the reason for the 
opens difference in the required bolt area is attributable to the neglect of — 
forces, some 33%, in the discusser’s method. ‘Nevertheless, there are 
_ significant differences in the bolt characteristics and design loads specified by " 
the national codes, which need to be examined,ie.; 42422 


A325(USA) 44.0(1. 70) 93.02.11) 
HSFG (UK) 35. 85.92.43) 43) 


' Bayes values are the equivalent stresses (in kips per square inch) based on 
the nominal shank diameter with the figures in brackets being the ratio of the 
at appropriate stress to the allowable stress. Note the ultimate /allowable stress 
_Tatio for HSFG bolts is 2.43 and not about 2.00 as suggested by the discusser. 
‘ If one now considers the required area for both types of bolt for a factored 


be of say 200 kips then using A325 bolts the required ann would need to 
be 2.67 sq in. as opposed t to 3. tee in. for HSFG bolts, giving a difference ; 


is s mainly due to the different aren stresses adopted for the bolts producing 

only some 8% difference in area between the two methods. The 33% owe 
force advocated by the writers is absorbed by the inherent reserve of rors 

& between proof and ultimate stresses of the HSFG bolt leaving an ae 


_ margin of safety against instant bolt fracture at plastic | limit load level. aa 
The discusser also comments that his method gives end-plate ices 
_ which are approx 90% of the values obtained from the writers’ recommendations. _ 


a Ref. 28 the discusser compares the end-plate thicknesses, determined for 


: different design criteria, for 36-ksi steel, knowing the A ,/A ,, ratio of the a 
_ sections. In the range of universal beam sections, equivalent tc to WF sections, 


‘method virtually gives the s same values as the split- -tee method. 
_ direct comparison can be made between the split-tee method and the tee-stub 
model used by the writers. The difference in the end-plate thicknesses can 
be accounted for by the definition of the effective span p,. The discusser achieves 
a smaller value of p, by assuming the application of the bolt load to be D/4 _— = 
_ nearer to the beam web, he justifies this as a result of extensive finite analyses 
(29). The location of the bolt load will always be a debatable point. Some 
United Kingdom commercial interests advocate that the value of D/4 should 
be (3/4) D or even D because it results in thinner end-plates and presumably - 
lower costs. Unfortunately such thinner end-plates carry the penalty of increased | 
im deformations within the connection and the subsequent loss of moment capacity a 
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DISCUSSION 


as stated earlier. Packer and Morris (14) illustrated that the local influence 


of the bolt head can be misleading when considering the general flexural behavior y 


of a plate element. The difference in end-plate thicknesses approximates to 
_ about a 5% drop in the ‘“‘moment”’ capacity of the bolts. 


_ By using smaller bolts and thinner end-plates there would be a tendency ; 
‘for increased connection deformations causing relaxation on the tension side 
of the beam and a lowering of the moment ‘capacity of the bolt group. This — 


m _ than those proposed by the modified split-tee methods are used in tests, then — 


_ the moment capacity reduces significantly. Such a reduction in moment capacity 
for a decrease in end-plate thickness was illustrated by the writers in a seocemt 
paper (33) in connection with lack of fit in flush end-plate 

‘Incidentally, the test connections by the discusser | have smaller cross” 
“centers of holes than would be the normal practical dimension used in the 
‘United Kingdom (about 5D). This automatically results in thinner end- -plates — 


; 4 becoming acceptable because the cross bending effect in the end-plates is reduced. 


When the two methods are examined in detail for the range of beam section 


“indicated, there appears to be only about 10% difference in both bolt area 
and end-plate thickness. This” really reflects. the objectives of the individual 


_ is confirmed by the discusser who states that when even thinner end-plates 5 


‘research programs, i.e., the discusser is interested in the elastic capacity of » 


a beam while the writers have concentrated on the er capacity with adequate 


_ As regards the discusser’s last two paragraphs concerning current work the 


_ writers cannot comment on the rather nebulous description of the work, except 


— to ‘Say that it den. be unrealistic to | assume a a fully active yield line extending — 
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32. _Agerskov, H., * 


Mann, A. P., and Morris, L. J., ‘““Lack of Fit in End- Plate 


the April 14-18, 1980, ASCE Convention and at Portland, 
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Dome EFFECT In Ri RC Sass: Ricw-Piastic At ANALYSIS" 


‘The writer read the paper 


‘Practical terms the central issue is the effect of unloading (or strain reduction) 7 
4 ‘on the stress block in the concrete compression z zone. : 


mess 


2a 


G. 7.—Typical Stress-Strain Relationship for Concrete, 
= T is of c course true that, when the position of the neu a axis moves towards 
_ the extreme compressed fiber material near the neutral axis undergoes a reduction 
ix of strain. The same thing occurs in an ordinary under-reinforced beam once 
_ the steel yields and, in this case, seems not to worry anyone. What is not clear 
is whether the strain reduction has any significant effect on the compressive 
force and it seems that this question is confused by the application of rigid-plastic 
7 shows typical stress-strain relationship for concrete and includes 
‘two unloading lines: one from a stage before the peak stress and one from Z 


NWI SLS, LS, England. 


Sune, 1980, by Mikael W. Braestrup (Proc. Paper 550). @ | 
be sehnic of Central London, 3) Maryieoone — | 


ZB a weapiealias zone » whose depth is halved while the peak strain doubles 
in a regime in which the product ‘‘peak strain x depth’? remains constant. 
UTRAL AXIS / F 
STRAIN / AXIS DEPTH 


8.—Development of S of Strains in Compression Zone alt 


STRESS 


~ 


STRAIN REDUCTIONS 


WITHOUT REGARD 


Finally Fig. 9 shows the stress Setstbations: for this compression zone at : 

the final stage as derived on the basis of Figs. 7 and 8 and also as ‘obtained 
directly from the envelope of Fig. 7 without regard to strain reductions. The 

differences between the two lines can be seen to be small. The reason is that 


the area where there is any appreciable strain decrement is one in which the . 
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— FIG, 9.—Final Stress Distribution in Compression Zone 
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_ the essential features of arching action rather better then a rigid-plastic — 
according to which yield of the concrete occurs over half the depth of the 


“of the load/ deformation properties of concrete in the rising branch. 


Dyn NAMIC 1c REDUCTION © OF STRUCTURAL Mope.s” 


Discussion by by Pip L. Gould,’ ‘Gould, ASCE 

The at author treats an subject which is of considerable practical 
_ importance. However, the development is incomplete since a firm theoretical 
basis for quantitatively assessing the effect of dynamic reduction has been 
established by Geradin (5) and this has apparently been overlooked. The latter 
_ reference incorporates the authors’ equations, for the most part, and goes far 
beyond in in providing a basis for error analysis. 
_ Some quantitative applications of Geradin’s technique to the finite element — 
wae of shells have been carried out by the writer and his co-workers. it 
was shown that the use of high-order approximations can detrimentally effect 
the calculation of eigenvalues: when kinematic condensation is igleyes, 


* Geradin, 4 M., ° “Error Bounds for Eigenvalue Analysis by Elimination of Variables,’ } 
Journal of Sound and Vibration, Vol. 19, 1971, pp. 11-132. 
«& Sen, S. K. and Gould, P. L., ‘Free Vibration of Shells of Revolution “a FEM,” 
Journal of the Engineering Mechanics Division, ASCE, 100, No. 
Paper 9629, Apr., 1974, 283-3038. 


DisTRIBUTIONS AND ExtreME Winp SpEeEps” 


‘The following corrections , should be made in the original paper; 


‘Page 2366, Eq. 3: should read [—In (1 — F)] instead of [—In (1 — os F) 
_ Page 2371, Table 2: the heading ‘‘Wind Speed” of Col. 7 and the eating 


October, 1980, by Charles A. Miller (Proc. Paper 15775), 


< 


-? Prof. and Chmn., Dept. of Civ. Engrg., School of Engrg. and Applied vant 
Washington Univ., Campus Box 1130, St. Louis, Mo. 63130. a Si. zane sate, 
a "December, 1980, by Emil Simiu and James J. Filliben (Proc. Paper 15909). fad ae'wDe 


member without there any strain OF Geliection. ihe writer thin hat 
| 
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“Data Sets” of Col. 8 should be diteted: Cols. 7 and 8 should have the common 
heading “Wind Speed DataSets.” 
Page 2371, Table 2, Cols. 4, 6, and 8: the numbers in the first, third, fifth, 
~ and seventh lines should be deleted; the numbers in the second and sixth lines — 
are percentages of sets best fitted by Weibull (y = 1) or Extreme Value Type 
y distributions; the numbers in the fourth and eighth lines are percentages 0 of 


sets ™ fitted by Weibull fy = = 2) or Weibull (y > 2) distributions. : one 
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